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Abstract 

The current theory of strong interactions, the quantum chromodynamics (QCD), is a 
non-abehan gauge theory, based on the gauge group SU(3). Despite its formal similarity 
to QED there are significant differences. It was shown in 1973 already that the coupling 
constant g increases for large distances pT|. This gave hope for the possibility to explain 



quark and gluon confinement. Soon after, in 1975, non-trivial solutions of the Euclidian 



Yang Mills equations were found, nowadays called BPST instantons [1^, which signifi- 
cantly influence the low energy structure of QCD. Many exact results are known for the 
1-instanton vacuum [|ll|, |18|, whereby an interesting phenomenological result of it is the 



explicit breaking of the axial ^(1) symmetry |]T^. On the other hand, a one instanton 

approximation, similar to a tree approximation in perturbation theory, cannot describe 

boundstates or spontaneous symmetry breaking. The next step was the analysis of exact 

T6| and approximate Jl^ multi-instanton solutions. There are two useful visualizations 



for these solutions. In one of these, instantons are interpreted as tunneling processes 
between different vacua. In the other interpretation, a solution describes an ensemble of 
extended (pseudo) particles in 4 dimensions. 

In conventional perturbation theory one computes fluctuations around the trivial zero 
solution. The correct quantization process is to consider all classical solutions of the 
field equations and their fluctuations. In the path integral representation of QCD the 
partition function is, hence, dominated by an ensemble of extended particles (instantons) 
in 4 dimensions at temperature g^. In the simplest case the partition function describes 
a diluted ideal gas of independent instantons. Unfortunately, this assumption leads to 
an infinite instanton density caused by large instantons, which obviously contradicts the 
assumption of a diluted gas. This problem is known as the infrared problem. The problem 



is avoided by assuming a repulsive interaction [^ which prevents the collapse. This is 
the model of a 4 dimensional liquid. Under certain circumstances the interaction can be 
replaced by an effective density. The Instanton Liquid Model in a narrow sense describes 
the QCD vacuum as a sum of independent instantons with radius p = (600MeV)~^ and 
effective density n = (200MeV)^. The correctness of this model is still being intensively 
investigated. So far the model is essentially justified by its phenomenological success. 

Numerical simulations of the Instanton Liquid Model allowed to determine a number of 
hadronic quantities, especially meson masses, baryon masses, hadron wave functions, and 
condensates [^ p^ . 

For computing the quark propagators and the meson correlators there are also analytical 
methods. The most important predictions are probably the breaking of the chiral sym- 
metry (SBCS) in the axial triplet channel |J25| and the absence of Goldstone bosons in 
the axial singlet channel. 

The largest part of this thesis is devoted to extending the analytical methods and to 
evaluating the results in (semi) analytical form. 

The meson correlators (also called polarization functions) will be computed in the Instan- 
ton Liquid Model in zeromode and l/N^ approximation, whereby dynamic quark loops 
will be taken into account. A spectral fit allows the computation of the masses of the a, p, 



cj, ai and /i mesons in the chiral limit. A separate consideration also allows computation 
of the 7]' mass. The results coincide on a 10% level with the experimental values. Further- 
more, determining the axial form factors of the proton, which are related to the proton 
spin (problem), will be attempted. A gauge invariant gluon mass for small momentum 
will also be computed. 

The thesis ends with several predictions which do not rely on the Instanton Liquid Model. 
In the 1-instanton vacuum a gauge invariant quark propagator will be computed and 
compared to the regular and singular propagator. Rules for the choice of a suitable 
gauge, especially between regular or singular, will be developed. A finite relation between 
the quark condensate and the QCD scale A will be derived, whereby neither an infrared 
cutoff, nor a specific instanton model will be used. 



Contents 



1 Introduction 



1.1 Methods to Solve QCD 



1.2 Contents 



7 

7 
9 



Theory of the Instanton LiquiH] 11 

2.1 Separating Gaussian from non-Gaussian Degrees of Freedom^ 11 



2.2 Effective QCD Lagrangiau in a Background Field *| 13 



O The Seniiclassical Limit *| 15 



2.4 Instantons in QCD 16 

2.5 (ojuarks| 17 

2.6 The Instanton Liquid Model] 18 



Light Quark Propagator| 20 

j.l Perturbation 'I'heory in the Muhi-instanton Background "^ 20 

3.2 Exact Scattering Amphtude in the one Instanton Background^ 21 

3.3 Zero Mode Approximation] 22 



EX 



Effective Vertex in the Multi-instanton Background '^ 22 

23 



3.5 A Nice Cancellation 



3.6 Renormalization of the Instanton-Density '*\ 24 

3.7 Selfconsistency Equation for the Quark Propagator "^ 25 



3.8 Some Phenomenological Results] 26 



3.9 Large N^. expansion^ 27 



3.10 Summary 2J 



4 Four Point Functions 



4.1 Introduction 



30 

30 

4.2 Large N^ approximation '''] 31 

4.3 Sohition of the Bethe-Salpeter Equations '^] 33 



4.4 Triplet and Singlet Correlators *] 34 



4.5 Summary! 35 



5 Correlators of Light Mesons 



36 

5.1 Analytical Expressions] 36 

5.2 Analytical Results] 37 

5.3 Spectral Representation] 38 



CONTENTS 



3^ Plot & i-'it of Meson Correlators! 41 



The Axial Anomaly 

i'lie Mass o 



3.1 



JTlii 



R. 



43 

Mesonl 43 



JTI Measurement of the Axial i^brm i'actor^ 47 



3.3 Axial Singlet Currents fc Anomaly! 48 

3.4 'I'he Proton Spin and its Interpretation] 49 



3.5 Reduction of the Proton Form Factors to Vacuum Correlators 52 



3.6 'I'he Axial I'brm I'actors G"f/^(g)| 53 

|6.7 The Anomaly Form Factor A(q') *| 54 

nie Cluonic Form Factors KYfl,{q 



3.8 



56 



|6.9 Discussion] 57 



7 Gluon Mass 



59 

7. 1 Introduction] 59 

7.2 Gluon Propagator *] 60 



7.3 Propagator in Statistical Background '^] 62 

7A A Naive Estimate of the Gluon Mass *| 63 



7.5 Expansion in the Instanton Density ^ 64 



7.6 QCD Propagators *| 65 



7.7 Propagators for Small Momentum '*'] 67 

7^8 Zeromodes 1 68 

69 



7.9 Conclusions and further development^ 



S Gauge Invariant Quark Propagator 



70 

1 Generalities On the Choice of Gauge] 70 

j.2 A Natural Gauge] 71 



B.3 On the Gauge in Instanton Physics] 72 

g.4 'Fhe Quark Propagator in Axial Gauge '^] 72 



g.5 Effective Quark Mass^ 74 

g.6 The Quark Condensate] 76 



g.7 Summary] 77 



9 Conclusions) 



79 

j.l Summary 79 

9.2 Outlook! 80 



3.3 Acknowledgements] 80 



A Appendix 

KA Notations 



81 
81 

A. 2 Instantons in Singular, Regular and Axial Gauge 82 

A. 3 Averaging over the Instanton Parameter 77] 83 



A. 4 Numerical Evaluation of integrals 84 

A. 5 Numerical Evaluation of the Convolution! 85 



B Figures 



86 



CONTENTS 



B.l Panorama .bunctionl 86 

B.2 Constituent Quark Mass in Regular, Singular and Axial Gauge] 86 

B.3 Pseudoscalar Triplet Correlator (tt)] 86 



B.4 Pseudoscalar Singlet Correlator {ri')\ 86 

B.5 Scalar Triplet Correlator {d)\ 86 

p. 6 Scalar Singlet Correlator (a)] 86 

B.7 Axial Vector Correlator {ai, fi) 



86 



|B.8 Vector Correlator (p, u;) 86 



C References 



95 



The chapters marked with * are more technical in nature. 



List of Tables 



3.1 Dependence of various quantities on the parameters of the instanton liquid 



model nB.,p,NcAQB.)- I 28 

5.1 Chiral Symmetry Breaking and Superconductivity) 39 



5.2 Mesonic correlators 



40 



5.3 Meson mass m^, coupling constant A» and continuum threshold E^ obtained 



within the instanton liquid model in this work {1/Nc expansion), froni 



numerical simulation and from experiment. | 42 

|6.1 Proton form factors at zero momentum transfer) 52 

A.l Asymptotic behaviour of ^y^(p) ) 83 



Chapter 1 



Introduction 



The fundamental theories which describe the observed interactions are all gauge theories. 
They describe gravitation, electroweak and strong interaction. In the quantized version 
the forces between particles (fermions) are mediated by gauge bosons. QCD describes 
the strong interaction between quarks and gluons. It is a SU(3) gauge theory with the 
LagrangianQ 

C = ^G^.Gr + E ^^(^^ + ^^^)^■ (1-1) 

Taking into account only the light quarks u, d and sometimes s and setting their mass 
to zero, one arrives at a theory with only one parameter, the gauge coupling constant g, 
which actually is no parameter due to dimensional transmutation. Despite the fact that 
QCD looks so simple (e.g. compared to electroweak interaction) it is very hard to solve 
this theory due to nonperturbative effects. 



1.1 Methods to Solve QCD 



To clarify the role of instantons in QCD let me first give a list of the most important meth- 
ods used to tackle QCD, starting with very general methods applicable to any quantum 
field theory (QFT) and ending with more specific approaches: 



Axiomatic field theory: Wightman/Oswalder&Schrader have stated a set of 
Minkowskian/Euclidian axioms for vacuum correlators a general QFT should re- 
spect (analyticity, regularity, Lorentz invariance, locality, . . . ). It is clear that the 
theorems derived from these axioms have to be very general because no Lagrangian 
is used PI. 



Throughout the entire work the Euchdian formulation of QCD is used |lJ 
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Haag-Ruelle/LSZ Theory: S-Matrix elements are related to vacuum correlators. 
The S-Matrix is the central object for particle phenomenology, containing such im- 
portant information as cross sections, form factors, structure functions, .... Vacuum 
correlators are more suitable for theoretical studies because they avoid the need of 
explicitly constructing the Hilbert space, which is an extremely complicated task 
beyond perturbation theory. 

Quantization: One might think that quantization should not appear in a list of 
methods for solving QFT because it is the method to obtain and define a QFT. On 
the other hand, besides canonical quantization there are other ways of quantizing 
a theory. The most popular is the path integral quantization [Q. In textbooks for 
particle physics it is usually only used as an abbreviation to derive theorems more 
quickly. In general (beyond perturbation theory), different quantization methods 
lead to different physical and mathematical insights and different methods to solve 
the theory. Variants of the path integral quantization are the random walk quanti- 
zation used in lattice theories, and stochastic quantization. 

(Broken) Symmetry: Every degree of freedom like spin, flavor and color is the 
possible origin of (approximate) symmetries like SU{2Nf) or subgroups and SU{3) 
gauge invariance. Conserved currents and Ward identities can be obtained. In the 
case of light quarks, one further has approximate chiral symmetry leading to PCAC, 
axial Ward identities, current algebra theorems, soft pion physics,. . . . Furthermore, 
QCD with massless quarks possesses an anomalously broken scale invariance, which 
is the origin of the huge field of renormalization group techniques [Q . 

Perturbation theory: Due to asymptotic freedom the coupling constant g decreases 
for high energies and perturbation theory in g is applicable. QCD can thus be 
solved for processes which involve only momenta of, say, more than 1 GeV. The 
small distance behaviour of vacuum correlators is, thus, calculable (x < 0.2 fm). 

Operator Product Expansion (OPE): An improvement of perturbation theory is 
to separate the small distance physics from large distance effects. The former is 
contained in the so-called Wilson coefficients calculated perturbatively. The latter 
nonperturbative effects are contained in a few vacuum or hadron expectation values 
of local operators which have to be determined from phenomenology or uncertain 



assumptions [|T3|. Vacuum correlators can be obtained up to distances of a; < 
0.3... 0.5 fm. 

QCD Sum rules: QCD sum rules are widely used to determine hadron masses and 
couplings. The general method is to assume the existence of certain hadrons and 
take a resonance+continuum ansatz in the Minkowskian region for some correlator. 
The Euclidian correlator is calculated from theory (OPE, lattice, instantons). Via 
dispersion relations one may match both in some Euclidian window by fitting the 
hadron parameters which leads to a prediction for them. 

Effective Theories: One may construct effective Lagrangians containing mesons 
and/or baryons more or less motivated by QCD or history or other physical branches. 
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Their parameters are determined from experiment or QCD as far as possible. 

A variety of phenomena have been explained qualitatively and calculated quantitatively 
with the methods listed above, but the large distance behaviour of QCD is still unsolved. 
There are at least two problems belonging to this domain: chiral symmetry breaking and 
confinement. Up to now chiral symmetry breaking was assumed and the consequences 
such as Golstone bosons were discussed within this assumption. In OPE one takes the 
nonzero values of the quark and other condensates from experiment but has no possibility 
to predict them from theory. The quark condensate is the order parameter of chiral 
symmetry and a nonzero values indicates spontaneous breaking of this symmetry (SBCS). 
Confinement has also to be assumed. 

Approaches to solve these problems are: 

• Lattice QCD: In principle the method is very simple. The continuum is replaced by 
a fine lattice covering a large but finite volume. The path integral is thus replaced 
by a finite number of integrals evaluated numerically. All vacuum correlators can 
be obtained for arbitrary Euclidian distances. Confinement and SBCS have been 
shown and other hadronic parameters are obtained. In practice, lattice calculations 
are much less straightforward than this sketchy description might suggest |^. 

• Instantons: As in lattice QCD one evaluates the Euclidian path integrals but 
now in semiclassical approximation. In addition to the global minimum of the QCD 
action A'f^ = Q used by perturbation theory there are many other local minima called 
instantons which have to be taken into account. Inclusion of light quarks leads to 
an effective 2Nf quark vertex responsible for SBCS. Although confinement cannot 



be explained, a lot of hadron parameters can nevertheless be calculated p7| , ^ 
suggesting that confinement is not essential for the properties of hadrons. 

Often only a combination of the results from various approaches allows contact to phe- 
nomenology. 



1.2 Contents 

In this work I want to give a quantitative and systematic study of the implications of 
instantons to hadron properties. For an elementary introduction into the classical and 



semi-classical theory of solitons and instantons I recommend |11 



Chapter ^ is an introduction to the semiclassical evaluation of nontrivial integrals. After 
developing the method in the finite dimensional case the partition function of QCD is 
considered and the instanton liquid model is introduced. 

In Chapter |^ the propagator of a light quark is calculated. The approximations which 
have to be made are stated and discussed. It is shown that for one quark flavor the 
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1/Nc expansion is exact. The constituent quark masses and the quark condensates are 
calculated for u,d and s quarks. 

The same approximations are used in Chapter ^to calculate the 4 point functions. Special 
attention is payed to the singlet correlator where a chain of quark loops contributes and 
is not suppressed in the large- A''^ limit. Within one and the same approximation we get 
Goldstone bosons in the pseudoscalar triplet correlator but no massless singlet boson. 

In Chapter [^ meson correlators are discussed and plotted. Employing a spectral ansatz 
it is possible to extract various meson masses and couplings. They are compared to the 



values obtained from extensive numerical studies of the instanton liquid |28| and to the 
experimental values. 

In Chapter |^ the axial anomaly is examined. The most interesting quantities in the axial 
singlet channel are the mass of the r]' meson and the spin of the proton. A combination 
of the axial anomaly and the scale anomaly with ideas from instanton physics ?7i^/ can be 
determined with an accuracy of 10%. The focus is on the discussion of the proton spin 
and its computation in the Instanton Liquid Model. 

In Chapter |^ the ghost and gluon propagator in the Instanton Liquid Model will be 
computed. For small momentum the ghost mass is 340 MeV and the gluon mass is 480 
MeV. The masses are independent of the gauge parameter ^. 

Chapter \^ attempts to predict several quantities without relying on Instanton Liquid 
Model. In regular gauge it is possible to derive a finite relation between the quark con- 
densate and the QCD scale A. Since the results depends critically on the choice of the 
right gauge, the chapter starts with a detailed discussion on the choice of a gauge and 
with a calculation of a gauge independent propagator. 

A summary of the results, new in this thesis, can be found in Chapter |^, followed by open 
questions, which deserve further investigation. The Appendices contain a list of used 
notation and explicit expressions for one instanton and its Dirac zeromode in singular, 
regular, and axial gauge. Furthermore, methods to compute Fourier transforms and 
convolutions of Lorentz covariant functions are described. 

All Figures can be found at the end of this work before the Reference section (sorted with 
respect to topics). 

The various chapters can be read independently of each other. Sections marked with * 
are more technical in nature. 



Chapter 2 

Theory of the Instanton Liquid 



This chapter gives an introduction into the semiclassical evaluation of the QCD parti- 
tion function. For this, the fluctuations around the classical solutions of the Euclidian 
equations of motion are separated into (approximately) Gaussian and collective coordi- 
nates (Section p.2| ). In Section ^]3| the approximately Gaussian degrees of freedom will 
be integrated out. The classical solutions of the Yang Mills equations can be classified 
with respect to their topology N & Z and are called multi-instanton solutions. For an 
introduction into the mathematically beautiful theory of instantons the reader should 
consult 1^. An elementary introduction into the classical and semiclassical theory of 



solitons and instantons can be found in [|T^. The classical one- instanton solution and 



its corresponding semi-classical partition function will be presented in Section |2.4| . The 
multi-instanton solutions can be reduced to the one-instanton case, in the case of well- 
separated instantons. Section ^3| describes how quarks modify the partition function. 
In Section ^.6| the infrared problem will be discussed and the Instanton Liquid Model to 
"solve" this problem. This model is the basis of Chapters ^^ 



2.1 Separating Gaussian from non-Gaussian Degrees 
of Freedom * 



In the path integral formulation solving QCD or any other QFT is equivalent to calculating 
the partition function Z including external source terms. Of course this is not a simple 
task because one has an infinite number of degrees of freedom. One way to tackle this 
problem is to reduce the number of degrees of freedom by integrating the (nearly) Gaussian 
degrees perturbatively leaving the non-Gaussian degrees, called collective coordinates, to 
be treated with other methods. 

Let me first describe this method most generally without referring to QFT. Our aim is to 
calculate 

Z= /"d"a;e-^[^l (2.1) 

11 
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where S is some real valued function depending on the n dimensional vector x G M"'. Let 
the integral be dominated by values of x which lie in the vicinity of a fc dimensional sub- 
manifold of iR" which may be parameterized by x = /(7), 7 G M^. Vectors in the tangent 
space of /(7) are called (approximate) zeromodes, because S'[x(7)] is (nearly) constant in 
this direction. Usually /(7) represents some degenerate or approximate minima of 5*. In 
addition 7(7) may contain points which do not contribute much to the functional integral, 
this will cause no error, but f{'~f) must not forget any significant points. In figure [B.l 



a two dimensional example is shown containing a river (one dimensional manifold) with 
steep mountains aside his banks. So the integral will be dominated by the shaded area. 
This area can be parameterized in the following way: 

x = f{l) + y , x,|/GiR" , ^EM" . (2.2) 

To make this representation unique, we have to demand k linear extra conditions to y: 

y.g,(^)=Ci , l<l<k (2.3) 

E.g., Cj = and (74(7) = df{j)/d'^i would fix y to be orthogonal to the "river" or equally 
stated: It disallows fluctuations in the non-Gaussian zeromode direction. Now every 
point X (at least in the vicinity of the river) can uniquely be described by y satisfying the 
condition (|2.3| ) and by 7 via ( |2.2| ). All we have to do now is to represent Z in terms of the 
fluctuation vector y and the collective coordinates 7. A convenient way is to introduce a 
Faddeev-Popov unit 



l = yAA5'(y^^(7)-Q)5"(/(7) + 2/-a:)$(x) (2.4) 

which serves as a definition of $. Inserting this into (|2.1|) and integrating over x yields 

Z = jd'^dry5\yg,{^) - q) $(/(7) + y) e-^l^^+^J . (2.5) 

The ^/-integration in ( p.4| ) can trivially be performed: 



^-\x)= dS5\{x-f{i))g.{i)-c,) . (2.6) 



For X = f{'y) + y the (5-function in the integral only contributes for 7' = 7. Thus we can 
expand the 5-argument up to linear order in 7' — 7 

<^-\f{l)+y) = /A'5ME(?/-%S^-^^(7)-^)(7;-7,)) (2.7) 



\detij{y-djgi - gi-djf)\ 
with y-gi{'j) = Ci 



-1 



^In the following we will omit the absolute bars; when necessary they have to be reinstated. 
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Inserting $ into ( p.5| ) we get 



One may write Z in a slightly different form, usually used in QFT, because it is more 
suitable for semiclassical approximations. Z is independent of q and therefore, although 
the r.h.s. of ( |2.^ ) explicitly contains Cj, it is actually independent of it. We can smooth 
the (5-function by a further multiplication with 

l = (27rO"'^/' fce-^S'"-' . (2.9) 

The determinant can be written in the form 

detA= fd7]dr]e^^" (2.10) 



where r] are anticommuting Grassmann variables (ghosts). Inserting ( |2.9|) and ( p.lO| ) into 
(p.8|) and performing the c-integration we finally get 



s,f = Y^y^(T.9^{l)9nlijy , (2.11) 

Sfpg = J2vt{y-djgth)-gth)-djf{-f))r]j . 

For a globally bijective transformation this representation of the partition function Z is 
still exact, for a locally bijective transformation the representation is exact in every order 
perturabtion theory, especially in semiclassical approximation. At this point the (approx- 
imately) Gaussian degrees of freedom y, t] and f/ could be integrated out in semiclassical 
approximation. The partition function then reduces to an integral over the collective 
coordinates. In Section |2.3| this step will be performed directly for the QCD case. 



2.2 Effective QCD Lagrangian in a Background Field 

Now it is time to return to QCD 

Z = jDA,e~'-^^^^ , SyM[A] = jdx^^Gl^G^,^ . (2.12) 

The background configurations which (approximately) minimize Sym will be denoted by 
^4^(7). A general gauge field can be written in the form 

A^ = {A,{^) + B,f , ^elR^ (2.13) 
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where -B^ are fluctuations around the background and A^ is a gauge transformed field 

A^ = SA^S^ + iSd^S^ , S = e'^ e SU{N^) . (2.14) 

As in the finite dimensional case we have to make this representation unique by introducing 
extra conditions, 

D^iA)B^ix) = Cix) , (2.15) 

to fix the gauge of the fluctuating field and 

jd^xijl{x-^)B^{x)=Ci , l<i<k (2.16) 

to avoid fluctuations in (approximately) zero mode direction. The derivation of an effec- 
tive action similar to (|2.11| ) can now be performed in full analogy to the previous case 
with only some notational complication. There is the following correspondence: 



finite example : i x y 'ji Qi 

QCD : i,x A B -fi,n{x) # 



(2.17) 



The Faddeev-Popov unit has the form 



1 = jd^-iDnDB,5{D^{A)B^)5\ ji^'^B^d'^)5{{A^ + B^f - A^)^A^] . (2.18) 

The steps to get an expression for $ are now: 

Add primes to 7, f2 and 5, insert A = A + B, linearize the last 5- argument around 
B' = Bfj_, perform the functional B' integration and linearize the remaining S arguments 



around 7' = 7 and Q,' = 0. Omitting the details of this calculation one gets p5 



^-\A{^)+B) = ld''yDn'5\Xi)6iY) , (2.19) 

X^ = Jd'xY.I^^',il)^-^^B,yi^-^,)+^;D,{A + B)^^ , (2.20) 

Y = J2D,{A + B)-^{i^-^,)+D,{A)D,{A + B)n' . (2.21) 

j 0'^3 



From (|2.15| ), ( |2.16| ), ( p.20| ) and ( |2.21| ) one can read off the form of the partition function 
Z: 

Z = N{^) fd^-fDB^D7]Df]6\ /"^;5^rf1r)e-^«^°[^-^''''*'l , (2.22) 



2.3. THE SEMICLASSICAL LIMIT 



15 



Sqcd 
Sym 



Sym[A + B]- Sgf[A, B] + Sfpg[A, B, r], r]] 



d'x—G%G^^{A + B) 



S,f = ^jd'x{D,{A)B^ 



^9f 

Sfpg 



2e 



d'xijUj) 



dA 



(2.23) 



M 



9lj 






Vj 



dA,, 



+ j<fxn{x)D,,{A)D^{A + B)r]{x) . 

Sqcd does not depend on the gauge parameter Q. For this reason the Q integration can 
be absorbed in the normahzation factor N{C,). ri{x) are the usual ghost fields originating 
from the gauge fixing. For every extra condition ( |2.16 ) one gets an additional ghost 
variable rji. For A = and no extra condition {k = 0) the action given above just reduces 
to the usual QCD action including Faddeev-Popov ghosts in R^ gauge 



S,f = ^Jd'x{d,B,f , Sfpg = J 



d^x rj{x)dij,Df,{B)r]{x) 



(2.24) 



Note that the action ( |2.22|) is still exact with the non-harmonic degrees of freedom 7^ now 
separated from the hopefully more Gaussian ones, B^ and 1]. 

For small coupling g it is now possible to establish Feynman rules from ( p. 231) in analogy 
to the case with no background. For this one has to know the "free" gluon, ghost and 
quark propagator in a given background A. If A is a non constant field even this is a 
very complicated task in contrast to usual perturbation theory around A = 0. For a 
multi-instanton configuration explicit expressions for the gluon and ghost propagator are 



derived in fl8 . 



2.3 The Semiclassical Limit * 

Before developing perturbation theory to all orders it is wise to study the semiclassical 
limit where one keeps only terms up to quadratic order in the fields. In QCD (and many 
other field theories) this is equivalent to lowest order perturbation theory, but around a 
very nontrivial background! 

Up to now we have not specified ip^^. A natural choice would be ■j/'^ = dA^/d'ji to fix the 
fluctuations to be orthogonal to the zero modes. Somewhat more convenient is to bring 
-0^ in background gauge: 



€ 



'dA, 



with n such that £';.(A)V'J, = 



(2.25) 
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Furthermore we assume that A minimizes the gauge action Sym which is true for widely 
separated instantons thus neglecting linear terms Sqcd- 

Up to quadratic order in the fields one has 

Sqcd = Sym[A] + jd^x^^B,K^M)Bu + jd^xn{x)D\A)7]{x) 

+ Y.V^%1^^V3+ d'^T.^{x)D,{A)^^,+0{^e\<l') , (2.26) 



1. 



K^, = -DH^, + 2iG^, + {l--)D^D, , G^, = F'Gl, , (F% = z/,,, .(2.27) 



Performing the integration over gauge fields and ghosts one gets an effective action de- 
pending only on the collective coordinates 7j: 



Z = f d'^j e-^^"^"'^ (2.28) 



dAA Bet{-D\A)) 



e-^^//W = det., ^;(7)^ ) tSWtI^^""^"'^' (2-29) 



The ^-function in ( |2.22| ) causes a restriction of the gauge field fiuctuation to be orthogonal 
to ip'^^. K'^^ is defined as i^^j, projected to the space orthogonal to -j/^^, Det' takes into 
account all eigenvalues of K' except the k zeromodes caused by the projection. 



^u 



2.4 Instantons in QCD 



The solutions of the classical Yang-Mills equations of motion can be classified w.r.t. their 
topology N & Z and are called N instantons^ solutions [1^. The 1-instanton solution 
has the well known form 

Ai^ix) = Oj 7],^, ^^ _ ^^^^ ^^ _ ^^^^ ^ ^^ m smgular gauge 

A%{x) = Ofvt,2 \^^''\2'l\ in regular gauge 

7/ = {zi,Oi, pj,Qi) = (location, orientation, radius, topological charge) 

The instanton parameters 7/ refiect the symmetries of the Lagrangian (translation, ro- 
tation, scale invariance, parity). It was a great achievement of 't Hooft to compute the 
functional determinant in the 1-instanton background. To get a finite results Zj has to be 



^ To simplify notations we will treat instantons and anti-instantons on the same footing. Both will be 
called instantons and are distinguished by their topological charge Qi = ± if necessary. 
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normalized to the A'^ = case, regularized and renormalized. The final result of the very 



complicated calculation [T^ for the partition function in semiclassical approximation is 

-) 

- Zq J reg 



fd^iD{pi) = ^ E fd%dOjdpjD{pj) = vJdpD{p) = V,D 



D{p) = ^Sr^e-'^^'^ = Cr^^p-'S^^ipAf 

4_gg-1.679JV. 



7r2(iVe-l)!(iV,-2)! 

7r2 87r2 1 

— , Slip) = -TT^ = ^l^-T ' 
9o 9i{p) pA 



So = ^ , Slip) = -Tr\ = ^^^~r ' A = Apy 



iTX 



2 



1 h' , I „, I 



S2{p) = -TT^ = b\n— + -\n\n— + 0{—-^) , 
glip) P^ b pA Any 

D{p) is the density of instantons of radius p, g{p) the running coupling constant in 1/2- 
loop approximation, b = yA^c and b' = ^N^. Sq is the classical action of an instanton, 
and Qq is the bare (unrenormalized) tree-level coupling constant. One can obtain an 
estimate for qq by replacing g^ with the running coupling constant at a suitable scale. 
This unlucky situation can be avoided by using the 2-loop expression for -D(p), i.e. by 
replacing 5*0 with Si and 5*1 with 5*2. But this replacement is an improvement only for 
small coupling. When p reaches the QCD scale A one should confine oneself to the tree- 
level or 1-loop approximation in order to obtain useful results, since the perturbation 
series is only asymptotically convergent. 



2.5 Quarks 

Additional fields coupled in a gauge invariant way to the gluon field can simply be incor- 
porated by adding the appropriate Lagrangian with gauge field A replaced hy A + B and 
performing the functional integration over the new fields. So every quark contributes an 
extra factor 

fo^D^e-I'^''^^'^'"'^'^ = Bet{ip+ im) (2.30) 

to the partition function Z where 

zD^ = id^ + A^ + B^ (2.31) 

is the covariant derivative. In the semiclassical approximation i?^ can be set to zero. D{p) 
has to be multiplied by the fermionic factor 

^, , f 1.34mp(l + m^p^ln(mp) + . . .) for mp < 1 , ^. 

and b and b' are now 

6 = lliV^_2 5' = l!jv2 _ l^iV^AT. + 1^ . (2.33) 

33^ Z ^ Z ^ 2N, ^ ' 
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2.6 The Instanton Liquid Model 

The sum of well-separated instantons {Qi = ±1) is also an approximate solution of the 
YM equations: 

N 

A = Y.Aj , 5[A]^iV^o 

7=1 

The partition function of this so called instanton gas is 

OO 1 

N=0 ^^ ■ 

The sum is dominated by instanton configuration with density N/V = D. Unfortunately 
D is infinite and the assumption of a diluted gas turns out to be wrong. The probability 
of small size instantons is low because D{p) vanishes rapidly for small distances. On the 
other hand for large distances D{p) blows up and soon gets large. This is the origin of 
the infrared problem which made a lot of people no longer believing in instanton physics. 



Those who were not deterred by that have thought of the following outcome |2^. For 
larger and larger distances, the vacuum gets more and more filled with instantons of 
increasing size. At some scale the instanton gas approximation breaks down and one has 
to consider the interaction between instantons which might be repulsive to stabilize the 
medium. The stabilization might occur at distances at which a semiclassical treatment is 
still possible and at densities at which the various instantons are still well separated objects 
- say - not much deformed through their interaction. So there is a narrow region of allowed 
values for the instanton radius. This picture of the vacuum is called the instanton liquid 
model. The idea has been confirmed in the course of years by very different approaches: 



Infrarot- Cutoff |T4| 
Hardcore assumption 
Variational Approach ||2 
Numerical studies p7[ 



Phenomenological success |^ 



The simplest suggestion is to introduces a cutoff pc and to ignore large instantons ||14|| : 



D,^= rdpD{p) 

Jo 



The cutoff has to be chosen sufficiently small such that the space time fraction / filled 
with instantons is smaller than 1 in order to justify the model of a diluted gas. 

This simple cutoff procedure can be improved by introducing a scale invariant (hardcore) 
repulsion between instantons, which effectively suppresses large instantons ||2^. This 
procedure has the advantage of respecting the scaling Ward identities which are otherwise 
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violated by the simple cutoff ansatz In |]25| such an repulsion has been found leading 
to a phenomenologically welcomed packing fraction. Unfortunately this repulsion is an 



artifact of the sum-ansatz as has been shown by |20|. Therefore the infrared problem is 
still unsolved. 

Nevertheless it is possible to make successful prediction by simply assuming a certain in- 
stanton density and some average radius. It seems that the vacuum can be described by ef- 
fectively independent instantons of size p = 600 MeV~^ and mean distance Lq = 200 MeV. 
The integral instanton density is fixed by the experimentally known gluon condensate |5j 



n = N/V, = l/Ll = ^ < Gl^C:^ >= (200 MeV)^,. (2.34) 

In Chapter ^ it is shown that this relation is not modified by light quarks. The ratio L^/p 
is estimated in different works to be 

{Lo/p)theor. = 3.0 . . . 3.2 (2.35) 

The Instanton Liquid Model is therefore defined by the following assumption on D{pi): 

D{pi) = nS{pi - p) , n= (200 MeV)^ , p = 600 MeV"^ . (2.36) 

The model describes very successfully the physics of light hadrons. Extensive numerical 
studies can be found in ^7\, |28[ . 

In high energy processes with momentum transfer p of 1 — 10 GeV, D{p) is often mul- 
tiplied with function, which is sharply peaked around p ~ p~^. The integral over p is 
then dominated by small instantons and is infrared convergent. The results are, hence, 
independent of the infrared cutoff. No additional assumptions have to be made. A sim- 
ilar phenomenon allows in Chapter ^ to derive a relation between the quark condensate 
and the QCD scale A without model assumptions. All other chapters are based on the 
Instanton Liquid Model. 



Chapter 3 

Light Quark Propagator 



In this chapter the average quark propagator in the multiinstanton background will be 



calculated. The methods developed in |2^ will be extended by taking into account dy- 
namical quark loops. In the first section the instanton background is treated as a classical 
external perturbation, but the background field is not small (e.g. in the coupling g) and 
therefore we have to sum up all Feynman graphs. This is possible in the case of one quark 
flavor within the so-called zeromode approximation. This is possible due to an exact can- 
cellation of certain graphs and by renormalizing the instanton density. As a byproduct 
we answer positively the question whether it is allowed to identify the gluon condensate 
with the instanton density in the presence of dynamical quarks. The quark condensate 



and a constituent quark mass are extracted from the quark propagator in section |3]^. In 
the last section it is shown that the case of two or more quark flavors can be reduced to 
the one flavor case in the limit oi Nc -^ oo. The results obtained in the one flavor case 
are therefore still valid when making this further approximation. 



3.1 Perturbation Theory in the Multi-instanton 
Background * 

It is well known how to calculate correlators in the presence of an external classical gauge 
fleld at least as perturbation series in powers of the external fleld A'^(x). In the case of 
QCD (or more accurately in classical chromodynamics) within the instanton liquid model 
the external fleld is a sum of well separated scatterers A = Y^j Aj called instantons with 
a flxed radius p and distributed randomly and independently in Euclidian space. 

For a while we will restrict ourselves to the case of one quark flavor and ignore gluon 
loops. The Euclidian Feynman rules have the following form: 

= Or) , 



A, 



X Ai 

— *-l 4— 



jjl + im 

4i ■ (3.1) 
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3.2. EXACT SCATTERING AMPLITUDE IN THE ONE INSTANTON BACKGROUND *21 

In operator notation, 

{p\So\q) = -t^Hp - q) , {x\4j\y) = Mx)6{x - y) , (3.2) 



{2ttY ' 

graphs are simply alternating chains of 5*0 and Aj. To average a graph over the instanton 
parameters one has to perform the following integration for each instanton: 



/4 J z, Q^^^ 1/4 



For example 




i^j 



'vi 



Jdjjd^j{p\So4iSo4iSo4jSo\q)MSo4jSo4i) ■ (3.4) 



The origin of the factor A^^ is the summation over all pairs of different instantons (/, J) 
yielding a factor A^(A^ — 1) ^ A^^. The quark loop is the origin of the minus sign and of 
the functional trace "Tr". 



3.2 Exact Scattering Amplitude in the one Instanton 
Background * 

Perturbation theory is suitable to study scattering processes. To achieve chiral symmetry 
breaking or bound states one has to sum up infinite series of a subclass of graphs or solve 
Schwinger Dyson or Bethe Selpeter equations. 



The first thing we can do is to sum up successive scatterings at one instanton 

xAj 



Aj X Ai xAi 

Vi := 4i + 4jSo4i + 4iSo4iSo4i + ■■■ = So\Si - So)s,' (3.5) 

where 

Sj = iS,'-4jr' (3.6) 

is the quark propagator in the one instanton background. 
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3.3 Zero Mode Approximation 

The quark propagator in the one instanton background can be represented in terms of 
the eigenvalues Aj and eigenfunctions ipi of the Dirac operator 

, Xi + im 



There is one zero eigenvalue iTj>il)i = (t Appendix |A.2| ) which makes the propagator 
singular in the chiral limit: 



im 

In the so called zero mode approximation one replaces the non zero mode part by the free 
propagator 

im 



s,= \iiim + sr^ . (3i 



5,-5o-^^^^^ . (3.9) 



Although this approximation is good for large as well as for small momenta it may be 
bad for intermediate ones, but what is more important is the fact that it is a wild ap- 
proximation and so might violate general theorems like Ward identities. In contrast, all 
other approximations we make are of systematic nature respecting all known symmetries 
of QCD. 

• semiclassical approximation (systematic) 

• multi- instanton background ("systematic") 

• large Nc expansion (systematic) (see next section) 

• zero mode approximation (wild) 

Note that every choice of a background gauge field is "systematic" in the sense of re- 
specting the symmetries of QCD as long as the background satisfies these symmetries on 
average. 

In the following sections we will see that the advantage of the zeromode approximation is 
so great that we cannot disregard this simplification. 



3.4 Effective Vertex in the Multi-instanton Back- 
ground * 

In Section ^^ the exact scattering amplitude Vj of the 1-instanton vacuum has been 
computed. In this Section a formal expression for the exact scattering amplitude Mj at 
one instanton in a multi-instanton bath will be derived. 

Let us consider a quark line with two scatterings at Vj and insert in between a number 
of instantons which differ from / and from all other instantons occurring elsewhere in 



3.5. A NICE CANCELLATION 
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the graph. This enables us to average over these enclosed instantons independently from 
the rest of the graph. Summation over all possible insertions with at least one instanton 
just yields the exact quark propagator minus the free propagator. Remember that direct 
repeated scattering at Aj is already included in Vj. 



Let us define 



+ -^Vi^Vi)- 



+ 



-(S>^<S)^<©- 



+ 



-^ 



■^- 



S-Sn 



(3.10) 



Mj = Vi + Vi{S-So)Vi + Vi{S-So)Vi{S-So)Vi + . 
= Vi + Vi{S-So)Mi 

This equation can be solved for Mj with the following ansatz: 



vLL 



Inserting Mj and Vj into ( p.lO ) we get 



li = m + i{^IJi\SQ {S - So)Sq \^j) 



(3.11) 



^S,'\^|JJ){^J\S^' = —(1 + -{HSo\S - So)S,'\^j))So'm{HS,' (3.12) 



(3.13) 



The 1-instanton bath causes a replacement of the current mass m with an effective mass 
/i. 

3.5 A Nice Cancellation * 



It is possible to arrange the graphs in such a way that every M/ occurs only once. Consider 
a graph containing two scattering processes Mj at the same instanton. The interesting 
part of the graph has the following form 



P 



q 



ipi{p)i'\{q) 



IfX 



i^i{r)ij\{s) 



ifj, 



(3.14) 



a.p/ip are the color/Dirac indices at the quark leg with momentum p. In a physical process 
in addition to the graph containing the above subgraph there exists another graph with 
only two quark lines interchanged. 

°pr7/)rfT)W/)iLQn«s arr7/)WrW/)U<7n«9 

(3.15) 



^/(p)^l(s)l"^"'-r^/(r)V'l(g) 



tjj 



tjj, 
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As usual, the interchange of two quark hues causes a minus sign in the amphtude. In- 
specting the two expressions, we see that they coincide except for the sign, thus there 
exists a complete cancellation 



Whenever an Mj occurs twice or more than twice in a graph there exists another graph 
with opposite sign. Both contributions cancel each other and can be ignored. So a quark 
can scatter only once at every instanton. This can be seen in another way: Because of 
Fermi statistics every state can be occupied only once, and there is only one state for each 
quark in the zero mode approximation, namely the zeromode. 

There are two equivalent descriptions of Feynman graphs: 

1. Draw all topologically distinct graphs with non- numbered vertices and assign a 
symmetry factor to each graph, 

2. Draw all topologically distinct graphs with numbered vertices and assign a factor 
1/V, where V is the number of vertices. 

If all possible graphs containing Mj are allowed it is not difficult to see within the second 
description that they can really be paired as stated above. 



3.6 Renormalization of the Instanton-Density * 

Up to now the cancellation is incomplete because not all graphs are allowed. Consider 
e.g. 



+ -^-4^}—^ = 



not allowed not allowed 

As in the case of Vi both graphs are not allowed. Another example is 



L(iu^_J 



-«- 



not allowed most general tadpole 



3. 7. SELFCONSISTENCY EQUATION FOR THE QUARK PROPAGATOR * 25 

It is a general fact that all disallowed graphs can be paired with other disallowed graphs 
or with tadpoles and vice versa. 

What we have to do is to "disallow" all tadpole graphs. Every Mj can be surrounded by 
tadpoles which contribute with a universal multiplicative factor which can be absorbed in 
a redefinition of the instanton density ri/j. Using this renormalized density n^ the pairing 
is now perfect and the statement "every Mj occurs only once" becomes true. 

One further can show that in the presence of dynamical quarks this renormalized density 
has to be identified with the gluon condensate instead of the "bare" density because the 
same tadpoles contribute to the gluon condensate too. The often asked question, whether 
in the presence of light quarks the gluon condensate can be identified with the instanton 
density, can be answered, in the above sense, with yes! 

3.7 Selfconsistency Equation for the Quark Propaga- 
tor * 

Quark loops are no longer possible because they cannot be connected to another part of 
the graph via a common instanton. All graphs which can contribute to the propagator 
are chains of different MjS. 

— ^ — = ( — ' — + ^4?^ + ^4?^"*€^ + ^^@-*^-*^©-^/ 

\ / I^J^... 

The M'jS can be averaged independently 

M(p):=z(M,) = ^pV'(p) , (3.17) 

where ip' is defined in appendix 0. The resulting expression for the propagator now has 

the form 

M M M 11 

S = So + So— So + So— So— So + ... = {Sq^ + M)-^ (3.18) 

where M = M{p) is the momentum dependent mass defined above. There is just one 
thing to do: We have to solve the circular dependence 



nR -^m( 1 {M3| 



(Ei|) ^^ 



but yU is just a number, which makes the solution very simple. Inserting (|3.171) and (|3.18|) 
into (|3.13|) one gets the following equation for yU 



2ip'^{p)Mj 



/^ = ^ + /^ ,./(^';M)2 (p' + ^(^ + ^p)) (3-19) 



p2 + (m + Mp)2 
which may be solved numerically for different current masses m. 
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3.8 Some Phenomenological Results 

In the chiral limit ( 3.19| ) reduces to 

fj.' = nnjd^ ^^ = annp' + 0{nl) (3.20) 



p2 + Mp2 

/i^ is proportional to ra^ and thus M is proportional to ^/nR in contrast to a linear 
dependence on nn obtained from a naive density expansion. 

In the last expression the denominator has been expanded in the density and 

a = p-2 f d^p'^ip'^(p) = 6.6 (3.21) 

is a universal number. For the standard values of npt and p one gets 

pI = Q.Qurp'^ = (lOOMeV)^ , 
M{p = 0) = 7.7pv^ = 300MeV . (3.22) 

The exact solution of ( |3.20| ) which has been obtained numerically by iteration, differs 
from the leading density value by 15%: 

M(0) = 345MeV . (3.23) 



The momentum dependence of the quark mass is shown in figure p.2| . The mass m + M{p) 



may be interpreted as the mass of a constituent quark. At high energies it tends to the 
current mass, at low momentum chiral symmetry breaking occurs and the quark gets its 
constituent mass M(0). Note that this is not a pole mass but a virtual mass at zero 
momentum squared. 



Let us now take into account a small current mass m formally of the order y^nR. The 
selfconsistency equation now reads 



- + 4 + 0{nR) . (3.24) 



Solving it for p leads to 



A^o < Ai = ^"^ + \/t"^^ + /^o ^ m + po (3.25) 

For the strange quark p is increased by a factor of 2: 

p{ms = 150MeV) = 200MeV (3.26) 

It is interesting that rus + M(0) remains to be 300MeV. For zero momentum the increase 
of the current mass is just compensated by an equal decrease of the dynamical mass M(0). 
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From the propagator one can obtain the quark condensate 

(^^) := ^mtTcoiSix) - Soix)) = N, j d^iir>{S{p) - S,{p)) . (3.27) 

In leading order in the density one gets 

z(^^) = ^^^ = (G;:,Gn/32vrV • (3.28) 

/i 

This leads to the following condensates for u,d and s quarks: 

i{uu) = i{dd) = (250MeV)^ , (ss) = 0.5{uu) . (3.29) 

For heavy quarks there exists a similar relation 

i{^^) = {Gl^G^^) I ^%T^'^m + 0{m-^) , (3.30) 

which leads within 10% to the same value for the strange quark condensate. This nicely 
confirms the hypothesis that the strange quark can be treated as a light quark as well as 
a heavy quark. This hypothesis is used in heavy to light quark matching formulas. 



3.9 Large Nc expansion 



* 



Consider now the case of N ^ light quark flavors u,d.,s, . . .. The discussion of the one 
flavor case in the previous sections can be copied up to the pairing and cancelation of 
graphs which contain more than one M/ (|3.16| ). This is still true in the multiple flavor 
case but now both quark lines in ( p.l6| ) must have the same flavor because Mi always 



connects quarks of the same flavor. So we have the theorem: "every Mj occurs only once 
for each flavor" . From this point on the discussion of the one flavor case breaks down 
because there are now graphs contributing to the propagator containing quark loops. The 
simplest new contribution has the form 



'l^-M 



U ,C7\ U ^^ U 



Is this contribution small in some sense ? Yes it is ! Quark loops are suppressed by a 
factor l/Nc- In perturbative context this is extensively discussed in |^|, in instanton 



physics it was first used by [^. Although 1/3 is not a very small number the large N^ 
expansion seems to be a good approximation in various cases. 

Consider a graph and add to it a new quark loop consisting of 

N new instantons q ^ at 

S instanton scatterings ~ 
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Parameters 



Instanton density 
Instanton radius 
Number of colors 
Coupling constant 



Gluon condensate 
Quark condensate 
Constituent mass 
Quark mass 
Gluon mass 

Ghost mass 
Meson correlator 



Quark loop 
Instanton scattering 
Instanton occurrence 
Gluon loop 
Instanton scattering 



ur, p, Nc, {g 



R) 



n = N/Vi = urN, ^ (200MeV)^ 
p ^ (eOOMeV)-^ 

g = gn/Nc (not used in semicl. limit) 



{GO) = {G%Gnm^' 
(ip-^) ^ 0.39 N^p-^y/^ 
M{p) ~ py/nE 



MghosM 

(^r^(x)^r^(o)) 



10.9 pn](^ 
7.7 pn](^ ft 



=: riRNc 

« (253MeV)3 

SOOMeV 



420MeV 40 



trunc. 



300MeV|0 




1 -1/2 



NjlR 

iV2-l 
{Nl - l)-i 



Table 3.1: Dependence of various quantities on the parameters of the instanton liquid 
model nR,p,Nc,{gR). 



This multiplies the graph (see table |3.1| ) by a factor loop = 0{N^ 
cases are possible: 



1+N-S\ 



The following 



1 + N-S 



= 
<0 



<^=^ M = N <^==^ all instantons are new 

<^=^ the loop is disconnected 

<^=^ one old instanton <^==^ the loop is a tadpole 

the loop is suppressed by at least one l/Nc 



Disconnected graphs are cancelled by the denominator and tadpoles have been absorbed 
in ur. So quark loops are indeed suppressed in the limit iVc — *> oo. The same is true 
for gluon loops. This can be seen by the same argument using the Nc dependencies from 
table 13? 



3.10 Summary 



The quark propagator has been computed in the Instanton Liquid Model in zeromode ap- 
proximation and in l/N^ expansion. It was possible to absorb dynamical quark loops into 
a renormalized instanton density, which should be identified with the gluon condensate. 
In the case of one quark flavor the 1/Nc expansion is exact. The values for the dynamical 
quark masses of the up, down and strange quark have been computed. Dynamical quark 
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loops are suppressed in every order perturabtion theory in the hmit Nc ^ oo. For the 
quark propagator this fact remains vahd beyond perturbation theory. In the next chapter 
we will see that in meson correlators there are non- suppressed quark loops. 



Chapter 4 



Four Point Functions 



In this chapter the correlators for four quark fields in the Instanton Liquid Model in 
zeromode approximation will be computed. Of special interest is the observation that in 
leading order 1/Nc not all quark loops are suppressed. Summation of all graphs leads 
to Bethe Salpeter equations in Section W7^ which are solved in Section |0|. Section O 



summarizes the results for the connected and free flavor singlet and triplet correlators. 

4.1 Introduction 

In the last chapter we have derived comprehensive Feynman rules within the zero mode 
approximation: 



j/ + im 

A quark of a given flavor can scatter only once at a instanton / via Mj. Tadpole graphs 
are not allowed, they are absorbed in the renormalized instanton density hr which has to 
be used when averaging over the instantons. In the large Nc limit dynamical quark loops 
are suppressed and /i can be determined by ( |3.19| ). In the chiral limit 

/i2 = 6.6n^p2 + 0(4) . (4.1) 

In this section we want to calculate 4 point functions in the case of two quark flavors of 
equal mass in the limit N^^oo within the zero mode approximation: 



V- 

d{p- s + r - g)nrr'(P; s,q,r) = r 



"r' = (4.2) 



dxdydzdwe'^Py~'''+''"-'''\0\T^{x)T^{y)^P{z)T'^lj{w)\0) . 

The ip fields are u or d quark fields arbitrarily mixed. Without restriction to generality 
we have taken the correlator to be a color singlet . These 4 point functions can be used 
to study meson correlators (fchapter ^ and quark form factors (fchapter P) 
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4.2 Large Nc approximation * 

The most general graph for the quark propagator is a sequence of different instantons Mj, 
according to the rules above. Similarly the most general graph for the triplet correlator 
{{ud){du)) consists of two quark propagators each containing every insanton only once. 
However, the u and d propagator may contain common instantons, e.g. 



It is always assumed that the left and right hand sides of the graphs form color singlets. 
Non-common instantons can be averaged like in the propagator case to yield graphs of 
the form ^ ^ 



^.(Ah-^Id^S-^aX-i- 



(4.3) 
where the thick line represents the full propagator 

— ^ — = ( ' — + ^4^ + ^^^-*^®-^ + —^jr^^)-*@-^) 

\ I i^j^... 

It can be shown that non planar diagrams like 



t^y->^m4- 



are suppressed by 1/Nc, where again the left and right hand sides of the graphs have to 
form color singlets. So only the ladder diagrams shown in ( [4.3|) contribute to the triplet 
correlator. 

One might think that the mixed correlator {{uu){dd)) is zero because the graphs neces- 
sarily include quark loops, or that only the two loop graphs contribute — but this is not 
the case! Let us first state the result and then discuss it. Graphs contributing to the 
mixed correlator are chains of quark bubbles 

d 



(4.4) 
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Application of the Nc counting rules shows that this chain is of order l/Nc- Taking the 
color trace at the left and right hand side of the chain we see that the mixed correlator is 
of order Nc- Using ( |3.16|) it is clear that the triplet correlator is of the same order. 



What is wrong with the derivation of quark loop suppression in the last chapter ? The 
main assumption was that every graph containing a loop can be constructed from a graph 
not possessing this loop by simply adding the loop. Eliminating a loop from the bubble 
chain ( [4.4|) yields a disconnected graph, but we only consider connected 4 point functions. 
So the quark loop chain cannot be constructed in a way needed to prove quark loop 
suppression. 



In the case of meson correlators one can take another point of view. The disconnected 
two loop contribution is of order N^ but (except for the scalar case) the contribution is 
zero. So the bubble chain is a subleading graph of order A''^ and nothing has been said 
about the form of subleading graphs. 



Nevertheless, all connected graphs can be obtained starting from {\i.4\ ) by adding further 
instantons and bubbles — but now Nc counting rules tell us that every attempt results in 
a 1/A^c suppression. Therefore the bubble chain is the most general leading order graph. 



Nothing has to be changed for the correlator {{dd){dd)) except that chain (O) must start 
with d. 



To calculate the connected 4 point functions one must now average and sum up the chains. 
Alternatively this can be represented in recursive from usually called Bethe Salpeter 
equations: 
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(4.5) 
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4.3 Solution of the Bethe-Salpeter Equations 



* 



Before solving the BS equations we have to construct the kernel. The l.h.s. of the kernel 
always forms a color singlet because of the restriction to color singlet correlators. In the 
following we work in singular gauge suitable for small momentum (| Chapter ^. 

Contracting the color and Dirac indices on the l.h.s. and using the formulas of Appendices 
A.2| and |A.3| one gets 



P—f- 



-f—q 



r 



s —h^Ih-l— r 



-^^{^^i{r)4{s)^Tjli,i{p)4{q)4)i = (4.6) 



nR 
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The kernel can now be determined to be 
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TiRNc 



M,M,MrMJ{p - s + r - q){SX6X + lXll\Win. ■ (4-7) 



qi ULs ULq 



The result is just proportional to the nonlocal version of the 't Hooft vertex between color 
singlet states. 

The solutions of the BS equations have a very similar structure: 
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Aq and A^ are scalar functions depending only o\it=p — s = q — r. The proof is simple: 
The kernel has the structure ( [4.8|) with Aq = A^ = 1. The product of two vertices yields 
the same structure: 
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urN, 



MpMqMrMsSij) -s + r-q) 



{AoFoBoit)6l6l + A,F,B,{tHrjlrJ6Z6 



s Oiq 



(4.9) 






idpds) — M^MsiiD{S{p)Sis)) , 
ur 

{dpds)—MpMstiD{S{p)-i^S{s)-i5) 
riR 



(4.10) 



(dpds) := / dpdsd{p — s — t) 



In other words, the vertices of structure ([4.8|) build a closed algebra. The reason for 
this simple result is that the kernel is a simple product function up to the momentum 
conserving 6. 



Using ( [4.7| ), ( [4.8|) and ( ^4.9|) the BS equations ( [4.5| ) reduce to primitive algebraic equations 
for Go/5it), Ho/5{t) and Ko/^^t): 

Go/^it) = 1 + Fo/5(t)Go/5(^) , 
7^0/5 (t) = -Fo/5it)Ho/,{t) , 



with the solution 



G 



1 



H 



1 



K 



1-F ' " 1-F2 ' 1-F2 

where we have suppressed the index 0/5 and the argument t. 



(4.11) 



4.4 Triplet and Singlet Correlators * 



Because of isospin symmetry SU{2)f, mesons form triplets and singlets. Replacing tptp 
in ( [4.2|) by the triplet and singlet combinations (borrowing the notation from the pseu- 
doscalar correlator) 



71^ = —j={uu — dd)^ TT^ = ud^ -K = du, r] = —={uu + dd) 
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Therefore C^ = K — H = jz^- This coincides with G = j^ for the charged triplet 
correlator ((7r^)(7r^)) as it should be. In the singlet case we get C^ = K + H = —j^- 

Generalization to Nf flavors causes no problems. When adding propagators in ( [4.8|) to 
the external legs the final result for the connected 4 point function is 

A^ I 

n^°r""(P,s,g,r) = ''-JM,M,MrM,[C^{t)iiD{S{s)TS{p))iin{S{q)TS{r)) + 

nR " 

C5it)trDiSis)TSiph,)tTDiS{q)TSirh,)] (4.14) 

Nf-1 
Co/5(^) = ~ fiv - 1)F (t) + 1 ^°^ *^^ singlet correlator (4.15) 

CiiAt) = — for the triplet correlator 

°/'^^ Fo/5(t)±l 



Fo/5(t) are defined in ( |4.10|) . The correlators of a singlet with a triplet current are zero 
as expected. For completeness we give the expressions for the disconnected parts of the 
4-point functions. The following graphs may contribute, depending on the fiavor structure 
of the correlator: 

P ^ q 

r r =N,tTn{TSip)T'{s))^ip-q)^ir-s) 

s ^ r 

(4.16) 

p — I I — q 

r /\ ^/ r =-N^tTD{TS{p))tTD(T'{q))^{p-s)^{q-r) 

Note that the second two loop term is of the order A^^. However, as discussed above, in 
most applications it drops out or yields an uninteresting constant or only the connected 
part is considered anyway. 

For the triplet and singlet case one gets 



d{p-s+q-r)U'^'^r{psqr) = N^tTD{TS{p)T'S{s))d{p-q)d{r-s) + ' 



for triplet 

2 • (4.16) for singlet 
(4.17) 



4.5 Summary 

4 point correlators in the Instanton Liquid Model in zeromode approximation and 1/Nc 
expansion have been computed. The 1/Nc is the substitute for density expansion, which 
breaks down in the presence of light quarks. In contrast to the original "perturbation 
theoretic" 1/Nc expansion [^ not all quark loops are suppressed. In the fiavor singlet 
channel a chain of quark loops survives the Nc ^ oo limit. The obtained expressions for 
the 4 point functions will be discussed in the following chapter. 



Chapter 5 

Correlators of Light Mesons 



Meson correlators, also called polarization functions, contain information about the meson 
spectrum of QCD. Poles at p^ = m? show the existence of mesons with mass m. With the 
4 point functions, calculated in the Instanton Liquid Model (jChapter^), we also possess 
the meson correlators analytically, besides integration (Section |5 .11) . In Section ^]l]we will 



discuss them briefly, whereby the pseudoscalar channel is of special interest. Due to chiral 
symmetry breaking there are massless Goldstone bosons in the triplet channel (but not 
in the singlet channel), because instantons explicitly break the U{1)a symmetry. Both 
phenomenons have analogies in superconductivity. In order to compute the meson masses 
we make an ansatz for the meson spectrum in Section |5.3| . Comparing the theoretical 
curve with this approach in Section |5.4| allows to determine the masses of the lightest 
mesons in the various channels. The evaluation of the integrals and the fit will be done 
numerically. 



5.1 Analytical Expressions 

In the last chapter we have calculated various quark 4 point functions (f4.2|) . The meson 
correlators or polarization functions are just local versions of these vertices and can be 
obtained by simply setting x = y and z = w. In momentum space the meson correlators 
have the form 



nrr'(i) = n^*''(t) + n'=°""(t) = ^(^i^ J^= (5.1 

{dpds) j{dqdr)Ii{p,s,q,r) = - 1^X6**^(0 |Tjr(x)jr'(0)| 0) , 

jr{xy" = ^{uTu{x)-dTd{x)) , jr'(O)^/* = -^(nr'«(0) - Jrd(O)) , 

(dpds) = / dpd^sd{p — s — t) , (dqdr) = / dqd^rd{q — r — t) , 

t=p — s = q — r . 
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From the explicit expressions of the 4 point functions obtained in the last chapter one can 
get, up to integration, analytical expressions for the meson correlators. The following list 
is a complete summary of all formulas needed to evaluate the meson correlators: 

nf^-(t) = N,l{dpds)tTD{rs{p)r's{s)) , 

n™r(t) = -iv,(Co(t)r°(t)r°,(t) + C5(t)r^(t)r^,(t)) , (5.2) 

, , 1 + for singlet correlator 

Fo/5(t) ±1 — for triplet correlator 

r°(t) = -L kdpds) JMpMstrDiSip)TSis)) , 



1 

Ur 



Tl{t) = -L kdpds) JM:pMsiiD{S{p)TS{s)^^) , 
^/nRJ ^ 

Fo{t) = — f{dpds)MpM,tTn{S{p)S{s)) , 

Ur J 

F^it) = — kdpds) MpM,tTD{S{ph,S{s)^5) , 

Ur J 

d 

pip'ip) = 2Trpz—[Io{z)Ko{z) - Ii{z)Ki{z)]^=pp/2 , 

/i = m+fd^ ^'^/^^]^l {p' + m{m + Mj,)) , (5.3) 

J p^ + [m + Mp)^ 

urN, = (200MeV)^ , p = (eOOMeV)-^ . 

5.2 Analytical Results 

Performing the Dirac traces leads to the following expressions: 

4 r^ ^ .MpM,{±{ps) - MpM,) 



4 /■ 
Fo/r>{t) = ndpds 



urJ' (p^ + M^){s^ + M^) 



4 r JMpMsi±ips)-MpM, 

—= lidpds)- = = — 

n^J^ "^ ' {p^ + M^){s^ + M^) 



^'/lii) = -^l(dpds) ' ,; ;:-^;;,; , /,/^ , (5.4) 






4^ . _ ^ jMpM,{MpS^ - M^p^) 



n^J (p2 + M2)(s2 + M2) 



Mp = m + Mp 



All other vertices F are zero. Consider the one instanton vertex ( [4.7|) (the kernel). It 
contributes only to the scalar and pseudoscalar correlator. From this observation one 
may have predicted that the connected part of all other channels is small because a 
contribution has to be a multi-instanton effect. Indeed, they are all zero as seen above 
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except for the axial correlator. Due to an extra factor M ~ y^^ in the numerator of F^g 
the connected part of the axial correlator is suppressed by OiriRj therefore it is small as 
expected and will be neglected in the following. 

Furthermore we will restrict ourself to the chiral limit, taking 171 = 0. Using the selfcon- 
sistency equation ( |5.3| ) one can see that 

F^it = 0) = 1 (5.5) 

is leading to a pole at t = in the pseudoscalar triplet correlator due to the F^it) — 1 
denominator in (|5.2|) . This is the massless Goldstone pion one expects in the chiral limit. 
A more extensive discussion can be found in [^ . On the other side in the singlet correlator 
the minus sign is replaced by a positive sign and there is no Goldstone boson in this case. 
Thus, the (two flavor) rj' meson is massive ! Unfortunately we cannot make any reliable 
prediction of the rj' mass because the kernal is very repulsive in this channel and no 
boundstate is formed. There have to be other attractive forces, e.g. confinement forces, to 
built an rj' boundstate. Similar things happen in the scalar triplet channel (compare figure 
B.5| and [B.(j|). But the most important thing is that there is no massless pseudoscalar 



singlet meson which is an important step towards discussing the f/(l)yi-problem. 

Both phenomenons have their analogon in the BCS-theory of superconductors, namely 
massless excitons and massive plasma oscillations. The phenomenons as well as the struc- 
ture of the calculation can be nearly 1:1 mapped onto, of course with differences in details. 
The 't Hooft vertex corresponds to the phonon-electron interaction. Further parallel can 
be read from table ^.1| . An introduction to the BCS-theory can be found in 



It is interesting to see that in leading order in the instanton density Fo[t) = —F^{t), which 
leads to a massless pole in the scalar singlet correlator. Numerically the cr-meson indeed 
turns out to be very light. Numerically a light scalar singlet meson cannot be excluded. 



5.3 Spectral Representation 

To extract phenomenological information from the meson correlators we make use of the 
spectral representation 

dxe'P^{0\Tjr{x)jr'mO) = J^ da^ D{a,x)p{a') (5.6) 

where 

Pip') = i2nfj:^iP-1nmJrmn){n\fA0m (5.7) 

n 

is the spectral density and 

D{m,x)= d^p^—^ = —^-^{mx)Ki{mx) (5.8) 

J p^ + m^ Att^x^ 

is the free propagator of mass m in coordinate representation. We have chosen the coordi- 
nate representation of U to be able to compare the plots directly with lattice calculations 
and with numerical studies of the instanton liquid [p8|]. 



5.3. SPECTRAL REPRESENTATION 
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Instantons in QCD 


BCS-Theory 


Phnomen 


Chiral Symm.-breaking 


Superconductivity 


Boundstates . . . 


gg-pairs=7r-mesons 


e~e~-pairs=Cooper-pairs 


bounded by attractive . . . 


instanton-induced 
't Hooft interaction 


phonon-electron- 
interaction 


\ \ 




in the . . . 


pseudo-scalar channel 


scalar channel 


Orderparameter 


Quark condensate 


Density of superconducting 
electrons (Cooper-pairs) (p) 


Spontaneous breaking of 


chiral symmetry 


gauge symmetry 


leads to . . . 


Goldstone bosons = 
massless pions 


Excitons = 
massless density fluct. 


Vacuum polarization 
diagrams make the 
massless excitations 
massive 




for the Coulomb interaction 
=^ massive plasma osz. 


■■■(Q)- 


in the flavor singlet case 
=^ massive rj' 



Table 5.1: Analogy between chiral symmetry breaking in QCD and the BCS theory of 
superconductors ^6l\] 
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Correlator 
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1=1 


1=0 


Pseudoscalar 


n5 = 


= {3535) 


hb 


TT 


i 


Scalar 


ni = 


= (iiji) 


Id 


6 


a 


Vector 


^,,- 


= (J/^Jm) 


7/. 


P 


uo 


Axialvector 


U^ -- 


= (J^J?) 


7m75 


ai 


h 



Table 5.2: Mesonic correlators 

The spectrum consists of mesonic resonances and the continuum contributions. If one is 
only interested in the properties of the first resonance one might approximate the rest of 
the spectrum by the perturbatively calculated continuum. 

One might think that the disconnected part only contributes to the continuum and the 
connected part will yield the boundstates. But this is not the case. On one hand, Bethe 
Salpeter equations have bound as well as continuum solutions. On the other hand consider 
a theory with weak attraction between particles of mass m. It is clear that there is only a 
cut above 2m and no boundstate pole in the free loop. However in the exact polarization 
function only a small portion of the continuum will be used to form a pole just below the 
threshold because the attraction is only weak. The Euclidian correlator will hardly be 
changed. Therefore, assuming weak attraction, we can already estimate the boundstate 
mass from the disconnected part. Of course in this example we need not calculate anything 
because we know that the bounstate mass is approximately 2m with errors of the order 
of the strength of the interaction. 

Assuming that all other forces neglected in QCD so far, especially perturbative corrections, 
are small and attractive in the vector and axialvector channel, we can obtain boundstate 
masses although in these channels up to our approximation there is no connected part. 
But things are less trivial than in the example above because the quarks do not posses a 
definite mass and we have to inspect the correlator to extract the meson masses. 



Let us start with the scalar and pseudoscalar correlator. 
m* is coupled to the current with strength 



The lowest resonance of mass 



A. 



(o|ji/5(o)b) 



(5.9) 



The rest of the spectrum is approximated by the continuum starting at the threshold Et^. 
* means vr, 7], 6 or a (see table p.2[ ). E^, is typically of the order 1.5 GeV and therefore 
the continuum can be calculated perturbatively. The spectrum thus has the form 



Pi/5(-s) = Xl^is -ml) + 



3s 



57r" 



Bis - El 



(5.10) 



Inserting p into (|5.6| ) one gets 



E{E„x) 



n'^x 



1/5^ 



n(v,(x) 



\lD{m^,x) +E{E^,x) 



3 {E.xf (27^3(^,0;) + {E,x)K,{E,x)) ^ 



16 



n'^x^ 



(5.11) 
(5.12) 
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In the next section m^, A* and E^, are obtained by fitting the phenomenological ansatz 
n'^**(x) to the theoretical curve 11'^'"™ (x) in the Euchdian region where the theoretical 
calculation is reliable. 

Consider now the vector and axial vector correlator. The vector current is conserved, thus 
the correlator is transverse and only the vector meson can contribute. In the chiral limit 
the same holds true for the axial current. In the singlet channel one has to be careful 
because there are two currents. A conserved one and a gauge invariant one which contains 
an anomaly. Up to now we have only calculated the correlator of the conserved current. 
Nevertheless to leading order in the instanton density the two correlators coincide and 
should be both conserved. 

For conserved vector and axial currents the spectral function is transverse: 

P^^Kp') = i-S,. + ^)p?ip') . (5.13) 

The coupling of the vector and axial meson to the current is given by 

tX.e, = {0\jI!\0)\p) (5.14) 

where e^ is the meson polarization. The spectral and polarization functions have the form 



3s 



pS(^) = ^Xl5{s - ml) + —e{s - Ef) , (5.15) 



Here * means p, u, ai or /i (see table p. 21 ). 



5.4 Plot & Fit of Meson Correlators 



The meson correlators are shown in figure |B.3| - p.8| , normalized to the correlators in 
lowest order perturbation theory: 

n?/5(-) = ^ ' "2?^ = -^ • (^-^^^ 

The diagrams therefore show the deviation from the perturbative behaviour. The nu- 
merical evaluation of the integrals are discussed in Appendices [A.4| and |A.5| . The meson 



parameters obtained by fitting the parameter ansatz to the theoretical curve are summa- 
rized in table |5.3| . Shuryak & Verbarshot |^ have obtained the same parameters from 



a numerical investigation of the instanton liquid model. Their values are also shown in 



table |5.3| . The parameters of the vector channel coincide extremely well with ours. For 
the IT and a meson there is a large discrepancy in the masses but this is not surprising: 
We are working in the chiral limit thus the pion mass has to be zero. A similar argument 
holds for the a meson as discussed above. The couplings fit very well. 

The discrepancy in the axial channels can have various origins which are under investiga- 
tion. Alternatively one may directly compare the graphs. They coincide very well even 
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Meson /^(J^^') 


m^MeV] v^[MeV] ^jMeV] source 


n l-(0-+) 


508±1 1276±33 1/N^ 
142±14 510±20 1360±100 simulation 
138 480 — experiment 


T]' 0+(0-+) 


^0 ? ? 1/iVe 

7^ ? ? simulation 
960 ? — experiment 


5 l-(0++) 


^0 ? ? 1/A^e 

7^ ? ? simulation 
970 ? — experiment 


a 0+(0++) 


433±3 506±3 1446±20 l/N^ 
543 500 1160 simulation 
? ? — experiment 


P 1+(1-) 


930±5 408±4 1455±33 l/N^ 
950±100 390±20 1500±100 simulation 
780 409±5 — experiment 


cu 0-(l — ) 


930±5 408±4 1455±33 l/N^ 
111 simulation 
780 390±5 — experiment 


ai 1-(1++) 


1350±200 370±30 1050±80 l/N^ 
1132±50 305±20 1100±50 simulation 
1260 400 — experiment 


/i o+(i++) 


1350±200 370±30 1050±80 l/N^ 
1210±50 293±20 1200±50 simulation 
1285 ? — experiment 



Table 5.3: Meson mass m^,, coupling constant A* and continuum threshold E^, obtained 
within the instanton liquid model in this work (1/A^c expansion), from numerical simula- 
tion and from experiment. 



in cases where a spectral fit does not work very well like in the 5 and rj' channel. The 
conclusion is that the terms neglected in our analytical treatment, but included in the 
numerical study |^, are small and usually give an correction less than 10%. These are 
contributions from nonzero modes and higher order corrections in l/Nf.. This is again an 
example for the surprisingly high accuracy of the l/N^. expansion. In the case of strange 
quarks the nonzero mode contributions will become more important. 

Finally, one should compare the numbers with experiment. As far as known, these num- 
bers are also listed in table |5.3|. A general discussion of the meson correlators and com- 



parison with experimental results can be found in [28 . 



Chapter 6 

The Axial Anomaly 



A variety of predictions concerning chiral symmetry breaking can be made within the 
instanton hquid model. Although the 't Hooft interaction |]l^ explicitly breaks the U{1) 



axial symmetry, instanton models are up to now not too successful in describing quanti- 
tatively the axial singlet channel. The most interesting quantities are the rj' mass and the 
spin of the proton. 



In Section ^]T] we compute m^' by combining various techniques. The result does not 
depend on the specifics of instanton physics. 



Sections |0| , |673| and ^^ are an introduction to the proton spin problem. In section ^]5] 
the proton form factors are reduced to vacuum correlators of 4 quark fields by assuming 
independent constituent quarks. The axial singlet quark and gluonic form factors are 
calculated in section |6l6| , [6?^ and |6]^ by using the propagator and 4 point functions of the 
instanton liquid model. Gauge (in) dependence is examined. A discussion of the results 
and a comparison with ||53[ is given in section ^ . 



In Sections 6!^6.5 we exceptionally use Minkowski- Not at ion. 



6.1 The Mass of the r]' Meson 

In many channels a direct calculation of the meson correlators in the instanton liquid 

model and a spectral fit lead to reasonable results for the masses of the lightest mesons 

2E| . This method even works in the axial triplet channel because the model correctly 

describes spontaneous breaking of chiral symmetry. In the axial singlet channel a strong 



repulsion prevents the formation of a meson |^, |22[ . The conclusion is, that there is no 
massless Goldstone boson in this channel, but the mass of the i]' remains undetermined. 
In this letter I want to calculate the mass of the rj' by combining quite different techniques. 
With the help of 

• current algebra theorems for the rj', 

• 1/Nc expansion, 
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• instanton model, 

• scale anomaly 

we are able to relate the rj' mass to the pion coupling constant fj^ and the physical gluon 
condensate. 

In leading order in l/N^ it is possible to relate the t]' mass to the B dependence of the 
topological susceptibility^ cPE/dQ"^ of QCD without quarks ||43[] : 



<' = jr[-^) ' ^j^ = J d'x {0\TQ{x)Qm0),onn (6.1) 

Q(x) = ^ti,GG(x) , Q = fd^Q(x) e Z 

4:71 J 

Q{x) is the topological charge density and Q the total charge. This formula is derived 
by arguing, that for large Nc the topological susceptibility is dominated by the r]' state, 
utilizing the axial anomaly ^T| and the relation /^ = /^/, which is exact for N^ -^ oo. 

The next step is to relate the topological susceptibility to the gluon condensate (0|A^(0)|0): 



N{x) = —tTcGG{x) , A^ = fd^xN{x) 



In instanton models the gluon field consists of instantons of charge Q = ±1. The exact 
N instanton solutions {Q = N) are selfdual {G^^ = G^^) and 

(o|rg(x)g(o)|o),o„„ = (o|Tiv(x)iv(o)|o)eo„n • (6.2) 



The exact anti-instanton solutions {Q = —N) are anti-selfdual {G^i, = —G^u) and (|6.2|) 
holds too, because the two minus signs cancel. Unfortunately these exact solutions are 
not the most important contributions to the partition function. 

The dominating configurations are instantons and anti-instantons in mixed combination. 
The simplest model is a dilute sum A = J2i ^i of instantons of mixed charge. G^j^u is 
then approximately selfdual near the instanton centers, approximately anti-selfdual near 
the anti-instanton centers and small far away from any instanton. In leading order in the 
instanton density we have 

Gf,^{x) = ±Gf,^{x) (6.3) 

where the sign now depends on x\ Let us define N^{x) in the following way: 

N{x) = N+{x) + N4x) , Q{x)= N+{x) - N^x) 

Nj^{x) is a sum of bumps near the centers of instantons and N-[x) near the centers of 
anti-instantons. (|6.2|) has to be replaced by the relation 



(0|TQ(x)Q(0)|0)eonn = {0\rN{x)N{0)\0),onn ' 4(0|riV+(a:)iV_ (0) |0),o„„, (6.4) 



1 {AB)conn^{AB)-(A){B) 
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Assuming independence of instantons and anti-instantons {{N^N^) = {N^){N^)) the 
equation reduces again to ( |6.2|) . We will see that this is a crucial assumption. 



The next ingredient is the scaling behaviour of QCD. Classical chromodynamics is scale 
invariant and the Noether theorem leads to a conserved scale current. In quantum theory 
the scale invariance is anomalously broken (like the axial singlet current). Ward identities 
can be derived, especially @] 

Jd'^X (0|TiV(x)iV(0)|0)eonn = ^(0|iV(0)|0) , b = yiV, (6.5) 

Therefore in a self(anti)dual background the topological susceptibility cPE/dQ'^ is pro- 
portional to the gluon condensate: 

= ^(O|iV(O)|O) (6.6) 

The relation is still valid, when there are statistically independent regions of selfduality 
and self ant iduality, as discussed above. It is in fact sufficient to assure independence 
of the total instanton/anti-instanton number N±. 1 have checked (|6.6| ) by using the 



theoretical one loop instanton density D(p) calculated in |T^. Only the p dependence 
D{p) ~ p~^{pA)'' is important. Due to the infrared divergence it is necessary to introduce 
an infrared cutoff, but one has to assure not to break scale invariance. A minimal change 
is to introduce two cutoffs f± in the total instanton/anti-instanton packing fraction. The 
packing fraction is the spacetime volume occupied by the instantons and is a dimensionless 
quantity. Scale invariance and independence of instantons and anti-instantons are ensured. 
The partition function Z is 

Z= Y: Zj,^Z],_ (6.7) 






jdp, . ..dp^^ D{p,) . ..D{p^J /± - ^EpH 



A lengthy, but quite standard calculation of statistical physics, leads to 

where c± = c±{f±, b) are constants independent of N±. Differentiation of In Z w.r.t. lnc± 
two times leads to ( |6.6| ). The result is independent of fj-. An attractive interaction would 
lower the susceptibility (compared to the density). This can be seen in the following 
way: In the extreme case of a very attractive interaction, all instantons will be bound 
to instanton- anti- instanton molecules, thus A^+ = A^_ and Q = 0. On the other hand a 
repulsive interaction would increase the susceptibility: 

(fE 4 

— — < -(0|A^(0)|0) for attractive // interaction (6.8) 

d'^E 4 - 

> -(0|A^(0)|0) for repulsive // interaction (6.9) 



dOS 
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Therefore the violation of ( |6.6| ) is a measure of the II interaction. 
All this should be compared to Dyakonov ||2^, where the relation 

d^E/dQ^ = {0\N{0)\0) (6.10) 



has been derived, which is similar to (|6.6| ). In this work a simple sum ansatz A = J2Ai 
has been made. This ansatz leads to a strong repulsion between close instantons, which 
is the origin of the missing factor 4/5. The result on its own and the comparison with 



( |6.9|) shows that some repulsive interaction is at work. Verbaarshot |^0[ has shown that 
this repulsion is an artifact of the simple sum ansatz. Using the much more accurate 
and elaborate streamline ansatz he showed, that the interaction strongly depends on the 
orientation and the average interaction is about 14 times smaller than those obtained in 



25| . Therefore the best thing one can do today is to assume no interaction at all and 
cut the packing fraction at some small value. There is also a more general argument that 
the relation derived in ||25[ must be wrong. The topological susceptibility is of 0{N^) 



whereas the gluon condensate is of 0{Nc). (|0|) is consistent with N^ -^ oo considerations 
only due to the presence of the 4/6 factor. 

Let us now continue with the calculation of m^/ . The physical gluon condensate in the 
presence of light quarks is 

(0|iV(0)|0)p;,j,, = (200MeV)^ . (6.11) 

Due to the presence of light quarks it is reduced by a factor a < 1. 

(0|iV(0)|0)p,,, = a(0|iV(0)|0)^ir'^'^ (6.12) 

Combining ([Ol), (|6.2|), ( |6.5| ) and ( |6.12| ), we get the final formula for the rj' mass is 

2 4:Nf 12 

One can see that rri^i ~ l/iVc because fl and ur and a are independent of Nf.. The 
largest uncertainty lies in the determination of a. Using again the instanton model, a is 
the determinant of the Dirac operator with the current quark masses replaced by effective 
masses: 

a= n l-34m;=^^p ^ 0.4 . . . 0.7 (6.14) 

i=u,d,s 

We have set the effective masses to the constituent masses 

ml^f = m"/^ = 300 . . . 350 MeV , mf -^ = 400 . . . 500 MeV (6.15) 

and p to the value of the instanton liquid model (p = 600 MeV~^). This estimate is 
consistent with the estimate of [^. Inserting Nf = 3 and f-^ = 132 MeV into ( |6.13|) we 



get 

m^/ = 884 ± 116MeV (6.16) 

which is in good agreement to the experimental value of 958 MeV. This result in turn 
confirms the assumption, that the interaction between selfdual and anti-selfdual regions 
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is small. The large uncertainty in a prevents more accurate statements, but ( |6.10| ) can 
definitely be excluded. 

Using nirj' as an input we can determine the gluon condensate in pure QCD 

(y 

-^{tr^GGy"'^'"'''''' = (246 MeV)^ (6.17) 

47r 

where we have again set Nf to 3. 

The factor 4/6 in (|6.6| ) is the essential term to get the correct Nc dependence of r/i^/ and 
agreement with the experimental mass. The discussion has shown, that the rj' channel can 
be an experimental device for testing the independence of selfdual and antiselfdual regions 
in QCD, which is an assumption in the simplest instanton models. It might turn out some 
day that the details of instanton models are wrong but the assumption of independent 
self (anti) dual regions remain valid. 



6.2 Measurement of the Axial Form Factors 

The forward matrix elements of the axial currents 

s^Aip = (ps|^7^75^|ps) , ^ = u,d,s 

can be interpreted as the quark spin content of the proton, in a sense defined more 
accurately in the following sections. Three independent linear combinations of Aip have 
been measured, thus allowing to extract their individual values. 



From the neutron /5-decay, using isospin invariance, one gets |BU|, E^ 



a3 = gA = Au-Ad = F + D = 1.254 ±0.06 . 
From the octet hyperon /3-decay, using SU{3)f symmetry, one gets [^ p4 



VSas = Au + Ad- 2As = 3F - D = 0.688 ± 0.0035 . 

From the spin dependent structure function gf of the proton, which has been measured 
by EMC [g and SMC ||, one can extract 

T„= f g{(x)dx = -Am + -Ad + -As + 0(aj = 0.142 ± 0.014 (6.18) 

Jo 9 9 9 

where we have given the world average value. 

Of special interest is the quarkspin sum, which can be extracted from the values given 
above, 

AS^^ = J-ao = Au + Ad + As = 0.27 ± 0.13 (6.19) 
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where the 0{as) corrections have been included. It deviates significantly from the naive 
quark model value ASg^ = 1- This deviation is the origin of the so called spin problem. 
Further the large polarization of strange quarks in the proton 

As = -0.1 ±0.05 

is counter intuitive because this indicates a large strange quark content of the proton. 

Much more could be said about proton spin phenomenology and the experiments. For an 
introduction and further references see ^, ^ ^7\, ^ . We will now give a more thorough 
definition and interpretation of AS*^^ and other quantities, which we want to calculate 
within the instanton model. 

6.3 Axial Singlet Currents & Anomaly 

It is well known that products of operators at the same spacetime point are very singular 
objects. In order to make the expressions well defined one has to regularize and renor- 
malize the operator products. An anomaly appears, if this procedure breaks a symmetry 
of the theory. The most important ones are the breakdown of the scale invariance and 
the breakdown of the axial symmetry |6^. In the following we are interested in the ax- 



ial anomaly [^. The operator product which has to be regularized is the axial singlet 
current, 

JA^) = Y. ^(a;)7M75g(a;) (6.20) 

q£{u,d,s,...} 

which seems to be local, gauge invariant and conserved]^ Unfortunately after regulariza- 
tion one of the three properties is unavoidably lost. Therefore we can define two different 
local currents, a conserved (c) one and a gauge invariant (GI) one. The third GI, con- 



served and non-local current is discussed in |^ in connection with the U{1) problem. We 



will suppress the summation over quark flavors and write ip for the quark fleld operator: 

J^^{x) = lim-?/'(x + e)7^75Pexp ( 2 / dz-A{z)\il>{x) 
JU^x) = limV'(a; + e)7^75V'(a;) (6.21) 

The difference between the two currents is described by the anomaly current K^: 

K,{x) = ^e^^^^tr^'^iG^'^-lA^An , J^,^ = J^^, + K, 
d'^KJx) = ^^^ti.GGix) = a{x) (6.22) 

271 

d^JU^x) = 2mJ5(x) , Js = #75^ . 

m is the current quark mass and Nf is the number of quark flavors. Note, that the splitting 
of J^5 in a conserved and an anomaly part is gauge dependent. There are attempts to 

^ We will use the term 'conserved' even for uriq ^ 0. Sometimes this current is called the symmetric 
current in the literature. 
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define both uniquely on physical grounds [^ . The intention is to define Jf^g as the naive 
parton model spin and K^^ as some gluonic contribution. The proton matrix elements of 
the various currents can be expressed in terms of real form factors Gi,Ki, AandJ: 



{V'S'\J',M\VS) = uAp')[l,l5G'i{q')-q,l,G^,{q')]uM 

{p's'\K^{0)\ps) = Us'ip')[-f^,'j5Ki{q^) -q^-f5K2iq^)\us{p) (6.23) 

{p's'\a{0)\ps) = 2MiA{q^)us'{p'h5Us{p) 
{p's'\J5{0)\ps) = iJ{q'^)us'ip'h5Us{p) 

M is the proton mass and q = p' — p. From (|6.22|) one can derive the following relations 
between the form factors: 

G^ = G^ + Ki , G2 = G^ + i^2 

2 2 

G\ - ^jGl = ^J , Ki- \-K2 = A (6.24) 

^ 2M ^ M ' 2M ^ ^ 

Gf-^Gf = Zj + A , 
^ 2M 2 M 

where all form factors are evaluated at q. The last equation relates only GI quantities. 
In the next section we show that the form factor Gf^ at zero momentum transfer can be 
connected with the proton spin. 



6.4 The Proton Spin and its Interpretation 

The stress tensor T^^^, is conserved {d^jT^'^ = 0) symmetric and GI and can be constructed 
from the Noether theorem. The angular momentum density tensor M^'^'' associated with 
Lorentz transformations can be expressed in terms of T^^: 

M""" = x^T^" - x^T^"" (6.25) 

M can be decomposed in spin and orbital contribution of quarks and gluons 

^,.p ^ M!;;p^, + M^;:;,^ + M^;'^, + M^:,;,^-^G\x''g^^-x^gn + d{- 



KZb = -i^l^{x''d^-xPd'')i^ , Mi;::;,^=-E^'^P^ij^„^,^ = -J^i (6.26) 



M^% = -G^'''{x''d''-xPd')A, , M^';f^i^ = G^'PA'' -G^^'AP 
The last two terms in M^'^'' do not contribute to the angular momentum operator 

f = -e'^'' jd\ M°^'=(x) . (6.27) 
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Taking the matrix element of Jz in a proton state, where the proton is ahgned in z- 
direction and at rest we get the spin of the proton 

AJ=^{ps\Jz\ps) = y'^'{ps\M'^\0)\ps) , AT = {p, s\p, s) = d\0) (6.28) 

The total spin of the proton is with no doubt 1/2 and we get the sum rule 

AJ = ALg + -ATF^ + ALg + Ag = - (6.29) 

where (ALg, ^AS*^^, ALg, Ag) are the (quark-orbital, quark-spin, gluon-orbital, gluon- 
spin) contribution to the proton spin, defined as matrix elements of the various parts of 
M given above. Therefore The GI axial current measures the quark spin contribution to 



the proton spin. The space integral in |6.27| can be cancelt with the state normalization 
and we get in covariant notation: 

s^ATP' = {ps\j'^^{<d)\ps) =^ ATP' = G'{\Q) (6.30) 

In the naive quark model the proton consists of three quarks at rest. There is no orbital 
and no gluonic contribution to the proton spin. This leads to the Ellis- Jaffe sum rule A J 
= }^ATP^ = 1/2. In the real world the identification of AS*^^ with the proton spin is not 
correct, because J^f^ measures the spin of the (nearly massless) current quarks whereas 
the proton consists of three massive (~ 300 MeV) constituent quarks. Further in a model 
of non-interacting constituent quarks the axial current which measures the constituent 
quark spin should be anomaly free because the anomaly is due to the interaction with 
gluons. Therefore the conserved current J^^ might be identified with the constituent 
quark spin operator. 

s.AE'^ = {ps\r^,mps) =^ aj:^ = gi{o) = i . (6.31) 



From (|6.24| ) we get 

AT.^^ = AJ:" + Ki{0) (6.32) 

which can now be interpreted in the following way: The spin of the constituent quarks 
AS'^ are formed by the spin of the current quarks AS'-^-'^ and a rest —Ki{0), which contains 
orbital and gluonic contributions. The origin of these contributions is not the motion and 
interaction of the constituent quarks inside the proton, because the constituent quarks 
are noninteracting and at rest in the naive quark model, but due to the formation of 
massive quarks from massless quarks. Therefore (|6.32|) may be discussed for an individual 
"constituent" quark. Further the gluonic configurations which are responsible for the 
generation of the quarkmass also determine the value of -ft'i(O). 

E.g. in a BAG model a massive quark is formed by confining a massless quark to a 
sphere. The spin of the massive constituent quark is the sum of the spin (^AS*^^) and 
the oribtal |(1 — AS*^^) contribution of the current quark. The BAG, which might be 
formed by nonperturbative gluonic configurations, is responsible for the mass generation 
and indirectly for the orbital contribution. From analytical and numerical calculations 
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we know, that in the BAG model the constituent spin is sphtted into 70% spin and 30% 
orbital contribution when starting with massless quarks. 

Whereas ( |6.32| ) is rigorously true, the interpretation of ^AS*^ as the spin of a constituent 
quark and its value ^ is questionable. One reason is, that an axial current which describes 
massive constituent quarks is by no means conserved in contradiction to J^g. 

There exists another relation between AT,'^^ and the form factor A at zero momentum 
transfer. Before deriving this relation we have to give a short discussion about the order 
of limits and massless poles. The following limits are taken: The spacetime volume goes 
to infinity (V4 -^ 00), because the universe is actually very large, the current quark masses 
go to zero (m ^ 0), because the up and down masses are very small and g — > 0, because 
we are interested in the forward matrix elements. In principle the results can depend on 
the order of the limits and therefore they have to be chosen consistent with the physical 
situation. This means, that if in the real world e.g. g <C m we first have to take q —>■ and 
then m —>■ 0. Actually we are interested in the forward matrix element {q = 0) and m y^ 
in the real world and the order of limits just stated applies. Through the cluster theorem 
connected correlators in coordinate space have to decay to zero when the separation of 
two arguments tends to infinity. Therefore there are no 5(g)-peaks in momentum space 
and the order of limits g ^ and V4 ^ 00 can be taken at will. Because m'^l^ ^ 1 we 
have to take first V4 — i> 00 and then m ^ 0. In statistical physics this is a well known fact, 
that a spontaneous breakdown of a symmetry only occurs, when there is a small explicit 
symmetry breaking term and the system volume tends to infinity. In the final end one may 
remove the symmetry breaking term. In QCD chiral symmetry is spontaneously broken 
(SBCS) and the small current quark mass is the explicit breaking of the chiral symmetry. 
Therefore it is mandatory first to take V4 — i> 00 and then m -^ pB|. Therefore we can 
use the following order of limits 

limllimLlim (...)]}. (6.33) 

This justifies the usage of the infinite volume formulation from the very beginning. 

In real QCD there are no massless particles (m^ 7^ 0). Therefore GI form factors have no 
massless poles especially 

q^Criq') ?^° (6.34) 

From (|6;2|), (|O0| ) and ( [Oil) we get 

TD 

AS^^ = -J(0) + A(0) (6.35) 

This relation is true whether there are Goldstone bosons in the axial singlet channel or 
not. Experimentally we know that the lightest particle in this channel is the rj' with a 
mass of 958 MeV much too large to be a Goldstone boson. Therefore J(0) remains finite 
in the chiral limit and we obtain 

AS^^ = A{0) for m^O (6.36) 

Assuming the non-existence of the axial singlet Goldstone boson from the very beginning 
the order of limits is of no importance in deriving (|6.36| ). Note, that ( |6.36| ) is only true, if 
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AS^ 


— 1^ r^ci m T 


Ki =^K2+ A 


AS^^-^G^^=AS^+/^i=^J+ A 


Nfm = M 


1 


= 0+1 


0=0+0 


1 - =1+0 = 1 + 


m = 


1 


= 1+0 


A-l= -1 + A 


A - = 1 +A-1= + A 


Instanton 


? 


= ? + 


= 1 +(-1) 


1 - = ? + = +(-1) 



Table 6.1: The proton form factors at zero momentum transfer q"^ = in the naive 
constituent quark model (Nfm = M), in chiral QCD (m=0) and in the instanton-liquid 
model (Instanton) . Experimentally A is 0.27. 

we take m = niu = nid = nig, although all three masses tend to zero. Otherwise additional 
nonsinglet currents on the r.h.s. of ( |6.35| ) would survive the chiral limit |]55[ . 



Combining (|6.24|) , (|6.34|) and ( |6.32|) we can conclude that 



2MAS' 



?^G'2(?^)k2=o = -g^^2(g^)|g2=o 



AE"^ is given by the pole residuum of Gg. Because G^-'^ has no massless pole AE'^ is also 
given by the pole of —K2. These massless poles are called ghost poles and they may truly 
appear, even if there are no physical massless particles, because Gg and K2 are gauge 



dependent objects. Note that all other form factors defined in ( |6.24| ) are GI and therefore 
free of massless poles. 



Table 6.4 summarizes the values for the form factors at zero momentum transfer for the 



following three cases: 

• the naive quark model of non- interacting constituent quarks of mass m = M/Nf, 

• chiral QCD and the identification of AE^ with the naive spin value 1, 

• the instanton liquid model. 

In the following sections we will calculate some of the form factors for a single constituent 
quark in the instanton liquid model. 



6.5 Reduction of the Proton Form Factors to Vac- 
uum Correlators 



In this section we will calculate some of the form factors defined above in the instanton 



liquid model. To apply the methods developed in [^ we relate the form factors to vacuum 
correlation functions 



■ir^s'ip') 



d^'xd^ze'P^- 



■ipz/ 



{p's'\B(fd)\ps) = 
% - M){-t^, - M){0\Tr]{x)B{0)r]{z)\0) 



(6.37) 



Usip) 



M is the proton mass and -B(O) is an arbitrary local operator. r){x) is a local operator with 
the quantum numbers of a proton e.g. a product of three quark fields in an appropriate 
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spin and flavor combination |Q. Assuming that 77 (x) tends to a free proton field operator 



for infinite times the proton states can be reduced and (|6.37|) is just an LSZ reduction 



formula for composite fields. For our purpose the following form is more suitable 
{p's'\BiO)\ps) = Z,n,,(p')[ lim S-\p')TB{p',p)S-\p)]u,ip) 



Tb{p',p) = ld^d^ze'P'''-'P'{0\T7]{x)B{0)f]{z)\0) (6.38) 



S{p) = fd^x e'P''{0\T7]{x)fj{0)\0) = /^'^ + conUnuum 



The advantage of this form is, that the explicit knowledge of the mass M is not needed. In 
Euclidian calculations like lattice-, instanton- and OPE-calculations it is always difficult 
to extract pole masses. 

This form can also be interpreted as a spectral representation of the 3 point function. 
Inserting two complete sets of states into the 3 point function and taking the limit p"^ = 
p'2 — » M^ to select the proton state one can directly attain (|6.38 ). 



If e.g. -B(O) is a quark current, the 3 point function is a product of 8 quark fields, which 
is too complicated to be evaluated in a multi-instanton background. Let us assume that 
the proton consists of three nearly independent quarks. Then the main nonperturbative 
properties of the proton come from the formation of constituent quarks out of current 
quarks. The forces which confine the constituent quarks in the proton are assumed to 
modify the properties of the proton only in a minor way, except that the proton is then 
stable. This assumption is justified by the success of the constituent quark model. The 
form factors of the proton are therefore the sum of the form factors of the constituent 
quarks. 77 has to be replaced by a single quark field ijj of fiavor up or down and M must 
be replaced by the constituent quark mass. In this case it is even more important to use 
( |6.38| ) because one does not expect a definite pole mass for the quark propagator. Looking 
at the quark propagator in the instanton liquid model we see, that the ^ term remains 
unrenormalized and therefore Z^ = 1. For a constant constituent mass this argument 
would be rigorously true. For a running mass it is plausible that Z^ is still approximately 
one. This fact is true in all models of chiral symmetry breaking I know. A conservative 
estimate is 

0.7 <Z,^<1 (6.39) 

In the following we will set Z^ = 1 remembering that this not an exact statement. The 
results for all form factors have to be multiplied with Z^. 

6.6 The Axial Form Factors G'f/2(g) 

The form factor of the current jr = ip'i^ip of a constituent quark can be reduced with the 
help of ( |6.38|) to a 4 point function 



trcolTjAp'^P)^'] = d'^xd^z e^P'^-*^^trcD[(0|r^(x)^(0)r^(0)^(^)|0)r] = (6.40) 
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Urr'iq-p,q-p',P,p') 

The polarization functions Hrr are calculated are defined and calculated in the instanton 
liquid model in p^ and other works. For F = 7^75 the connected part of the 4 point 
function is suppressed by 0{n^ ). In leading order in the instanton density only the 



disconnected part contributes and we get 

T^ciip',p) = S{p'h,^,Sip) (6.41) 

Inserting ( |6.41|) in (|6.38|) and comparison with (|6.23|) leads to 



{p's'\J%\0)\ps) = uAp'h,-f5Us{p) (6.42) 

GfV) = l , Gf(g^) = 

Note, that Ilrr' was calculated in singular gauge, but the connected part is suppressed 
in any gauge and the disconnected part only depends on the propagators, which cancel 
out anyway. The form factors GfL{(l) ^''^^ indeed gauge invariant. The result coincides 
with a model of free massive quarks. Further we see that the current is not conserved. 
Conservation demands q^G2 = MGi, which is clearly not satisfied by (|6.42|) . In the one 



instanton approximation one can work from the very beginning with the effective 't Hooft 



vertex [|T^] which explicitly breaks the U{1) symmetry and therefore contains the anomaly. 



The result for the GI form factors ( |6.42| ), although not consistent with the experimental 
value, is up to now at least theoretical consistent. 



6.7 The Anomaly Form Factor A{q) 



* 



We will now calculate the anomaly form factor A. Using again the reduction formula with 
insertion of the anomaly current -B(O) = a(0) we have to calculate the 3 point function 
Ta{p, s). In the instanton model the field operator a(0) is replaced by a classical field a^(0) 
where A = J2i ^i is a multi instanton configuration inserted in a. In a given background 
A the correlator can be written in the form 

{0\Ti;{x)aiO)iPiz)\0)A = aAiQ){0\Ti;{x)^{z)\0)A = aA{0)SA{x,z) (6.43) 

where Sa{x,z) is the quark propagator in the multi instanton background A. The r.h.s. 
has now to be averaged over the collective coordinates 7/ of all instantons. Without the 
factor 0^(0) this is just the averaged quark propagator calculated in |^. aA{y) is 2Nf 
times the topological charge density at spacetime point y. In the vicinity of an instanton of 
charge Qj = ±1 the charge density has a positive/negative bump and is small elsewhere. 
Therefore aA{y) is only nonzero when there is at least one instanton near y. Let us 
fix exactly one instanton in the vicinity oi y = 0. The orientation and charge of the 
remaining instantons can be averaged independently, but when averaging the locations 
zj the domain near y has to be avoided. The next step is to assume 2 instantons near y 
and so on. The relative error we make by neglecting these further contributions and by 
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forgetting about the restriction on zi are both of 0(71^). In leading order in the instanton 
density we can therefore fix one instanton near y = and take only this contribution to 
0^(0) into account. The remaining instantons can be averaged as in the pure propagator 
case and the diagrams which have to be summed and averaged are the same except for 
the fixing of one instanton /. The propagator consists of a chain of instanton scatterings 
Aj (J = l...A^). Repeated scattering at this vertex is allowed. There are two cases: 
The first case is that all instantons left to all occurrences of instanton / are different 
to all instantons right to all occurrences of instanton /. In leading order in l/N^ all 
instantons in the middle section from the first up to the last occurrence of Aj are different 
to the exterior instantons. The instantons on the left and on the right can be averaged 
independently leading to averaged multi-instanton propagators. Averaging the middle 
section, but fixing / leads to the effective vertex Mj. The free part of the correlator in 
momentum space is therefore 

Tt^%p,s) = {2NfQj{zj) ^-@^),= 



= -2tNjQ{p - s)^MpMsS{ph,S{s) (6.44) 

Qi{zi) is the charge density of one instanton of charge Qj = ±1 and Q{q) = \{qpY K2{qp) 



its Fourier transform^ for Qi = +1. For p^ = s^ = M^ the term J MpMs is just the onshell 
mass M. Inserting ( |6.44| ) into ( |6.38| ) and comparison with ( |6.23[ ) we get for the free part 
of the anomaly form factor: 

Af''\q) = -NfQ{q) , A^'^'iO) = -Nf (6.45) 

The second case is, that there are common instantons to the left and to the right of 
instanton /. The connected part of the correlator and the form factor are 



Tr''{p,s) = (2NfQj{zj 



p 



f^ 



C 



PP Q£ 



(6.46) 



= -4{Nf - l)tQ{p - s)Q{p - s)F,{p - s)^MpM,S{ph,S{s) 

A^"--{q) = -2{Nf-l)Q{q)Cl{q)F,{q) , A™""(0) = iV; - 1 (6.47) 

For one flavor the connected part is zero as it should. For two flavors the result can easily 



derived by using the formulas of |^ . The total anomaly form factor for zero momentum 
transfer 

A(0) = A^'^^{Q) + v4'^°'^"(0) = -1 (6.48) 



^ I apologize for the overload of the symbol K: K2{qp) is a modified Bessel function, K^{x) is the 



anomaly current and ^^1/2(9) its form factors. 
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is independent of the number of flavors! This result is welcomed due to the following ar- 
gument: The form factors of the axial singlet currents j^s should not depend on any quark 
flavor which is not involved in the particle state. One expects that they are independent 



of Nf. Due to (|6.22| ) matrix elements of a{x) must then be independent of Nf too. But 
this is not obvious because a{x) is explicitly proportional to Nf and the gluonic field is 
not flavor sensitive. The calculation given above shows how the quark interaction cancels 
the free part, which is proportional to Nf, so that the total form factor is independent of 
Nf at least at zero momentum transfer. 



6.8 The Gluonic Form Factors Kfk{Q) 

Now we come to the calculation of i^i/2(0). The previous calculation can be copied with 
minor changes. a(0) has to be replaced by -ft'^(O). This in turn induces the replacement 

2Q{zi)^G,{zi):=^K^^^{Q) , 2Q{q)^G,{q) (6.49) 

IMf 

G^{zi) is -ft'^(O) where the gauge field is an instanton centered at zj of charge Qi = +1 
and Gf^{q) is its Fourier transform. In regular gauge we get 

G7i') = W^aMO) = - 'i'llXl (6.50) 

G;''{q) = -tq,p'K,{qp)'-^-2zq^/q' 



With this replacement in ( |6.44|) and (|6.46|) and comparison with (|6.23| ) 7^1/2(0) can be 
extracted: 

Kr{q)=0 , Jim|^i^r(g) = l (6.51) 

In singular gauge we get 

Gri^) = G7(^) + ^, , G7^(g)=G-^(g) + 2zg,/g2^0 (6.52) 

An apparent observation is, that the anomaly form factor K2{q) is gauge dependent and 
receives a massless pole in regular gauge. The reason for this is the gauge dependence of 
the anomaly current Kfj_ itself. One can show that the forward matrix elements i^i/2(0) 
are GI for small gauge transformations. A gauge transformation is called small, when it 
can be smoothly deformed into the unit transformation. On the other hand the gauge 
transformation, which transforms an instanton from regular gauge to one in singular gauge 
is large, because the regular solution can not be smoothly deformed into a singular one 
due to the singularity. 
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The next striking observation is that the relation 

Kiiq) - ^K^iq) = A(g) (6.53) 

is violated in singular gauge as can be seen from ( |6.52| ) 

Knq) - ^Kr + A{q) (6.54) 



Surface terms are the origin of this violation. For the derivation of (|6.53|) one has assumed 
the vanishing of surface terms. If one replaces the plane wave solution for the state by 
a wave packet, the state and therefore the matrix elements decrease sufficiently fast at 
spacial infinity and there are no surface terms. A experimental state is always a more or 
less localized wave packet rather than an exact plane wave. Therefore in regular gauge 
there are no surface terms and 

K?\q) - ^Kr = A{q) (6.55) 

is valid for all q. In order to work in singular gauge we have to choose a space-time 
manifold iR^\{0} to exclude the unphysical singularity. This small hole should not affect 
the physics at large distances. Therefore all coordinate space integrals are integrals over 
the domain iR^\i?£(0). Partial integration can now lead to surface terms at zero. The 
surface term is non-zero in the case of G**"^ as can be seen from (|6.52|) . This is the reason 
for the inequality ( |6.54| ). It is surprising that not the slowly decaying regular gauge field 



causes a surface term at infinity but the strong singularity at the instanton centers in 
singular gauge leads to surface terms and to a violation of ( |6.53| ). 

The following conclusions should be drawn: 

1. not to consider gauge dependent objects like Ki/2{q) at all or 



2. save the relation (|6.53|) by using regular gauge although this violates the philosophy 

of [ESl or 



3. modify relation ( |6.53[ ) by including the surface terms and be careful when performing 
partial integrations. 



In the following discussion we take position 2. 



6.9 Discussion 



Comparing the results for the form factors AS^^ = Gf^(O), G^^(O), A(0), Ki{Q) = 
Kl^^{0) and -^'2(0) = -ft'2^^(0) summarized in the last row of table ^]4| we clearly see that 
they are in contradiction. It is not possible to determine the remaining form factors in a 
way that they are consistent with ( |6.24|) and ( |6.34|) . The most obvious contradiction is 



AS 7^ ^(0). An opposite sign of the anomaly would at least be theoretical consistent 
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and would lead to the naive expectations. The only candidate for this violation of the 
axial ward identities is the neglect of the non-zeromodes. All other approximations respect 



the symmetries of QCD as discussed in Chapter 3.3 



Forte IQ has derived the relation AS + ^4(0) = in the instanton model in the case of 



one quark flavor in quenched approximation and density expansion. AS was identified 
with AS'^ and -^'2(0) was assumed to be zero (although not explicitly stated). Therefore 



A(0) = Ki{0) and from (|6.32|) one can arrive at the welcomed result AS = 



In section p.8| I have shown that the anomaly contributes to K2 and not to Ki. This is 
the first discrepancy. Further, in section ^]6| I have shown that AS has to be identified 
with AS*^^. This is the second discrepancy. It may turn out that the inclusion of non- 
zeromodes removes the discrepancies in a way, that leads to a phenomenological welcomed 
small AS*^^. In the one instanton approximation the inclusion is manageable and has been 



performed by P2| for the meson correlators. The consistent extension to the instanton 



liquid and to the quark form factors was not yet manageable. 

These problems might be compared to calculations of the rj' mass. A brute force method 
of calculating the axial singlet meson correlator and extracting m^/ by a spectral fit is 
not successful too. More elaborate arguments, given in Section |6]^ allowed a successful 
determination of m^/ . The instanton model was only used as a motivation for a selfdual 
model of QCD. Maybe the same model is able to solve the proton spin problem rather 
than a brute force calculation. 

It might also be possible that the spin problem can not be solved on the level of individual 
constituent quarks formed out of current quarks but is connected with a strong interaction 
in the axial singlet channel between different constituent quarks. This possibility in 



connection with instantons is discussed in 69 . 



Chapter 7 



Gluon Mass 



Up to now the instanton vacuum has been treated on the classical level. In this Chap- 
ter we compute the quantum fluctuations around the 1-instanton vacuum for the gluon 
propagator. 



Section \l.l\ gives a brief introduction to the problems and phenomenological conse- 
quences of massive gauge bosons. In the other Sections we will calculate the gluon 
propagator for different gauges. In background ^ = 1 gauge it has the simple form 



^au — ^ 9iiiy/iP ~ M{p) )■ ^^ Section 7^ we will extract the inverse propagator from 



the terms quadratic in the fluctuations around a background field. In Section |7]^ we will 
average this expression over relevant background fields and expand it for large momentum 
in a way to get the gluon mass M{p). In Section [7.4| we explicitly calculate M in the 



multi-instanton background. For large momentum, however, this mass is cancelled out 
by terms which have been neglected. Furthermore, the mass is gauge dependent. To get 
reliable results for large as well as for small momentum we make in Section [7.5| a cluster 
expansion in the instanton density. For this purpose we need the gluon propagator in the 
1-instanton background. The relevant formulas to construct this propagator are listed in 



Section 7.6. Because they are a bit lengthy we will restrict ourselves in Section [7.7|to the 



case of small momentum. 



7.1 Introduction 



The question, whether creation of a mass for gauge bosons is possible without gauge sym- 
metry breaking, is very old. The explicit introduction of a gluon mass into the Lagrangian 
causes several difficulties. To preserve gauge invariance the gauge fields have to be cou- 
pled to massless scalar particles I^D], ^, which decouple from physical matrix elements. 



The mass cannot be a constant, but must vanish for large momentum in order to ensure 
renormalizability of the theory (soft mass) . Both phenomenons (massless scalar and soft 
gluon mass) occur in every model which dynamically creates a gluon mass. It should be 
mentioned that the terminology of a mass does not necessarily stand for a pole mass, 
but for the selfenergy. A dynamical mass means a non- vanishing selfenergy at p^ = 0, 
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i.e. if the massless pole has disappeared from the propagator. In |Q a gauge invariant 
selfenergy has been defined and computed. Solving the Schwinger Dyson equations leads 
to a mass of 500±200MeK. Due to asymptotic freedom the mass vanishes logarithmically 
for large momenta Mgiuon{p) ~ (lnp^)~^^/^^. Nonperturbative arguments lead to a decay 
^(Inp^)^^/^^ W^ RH]. Finally it should be mentioned that a dynamical gluon mass is not 



in contradiction to confinement [p4| , ^ , ^ , [37|| . 



Apart from regularizing all infrared divergences the gluon mass has a series of phenomeno- 
logical consequences. A direct consequence are gluon balls. In the simplest picture in 
which glueballs consist of A^ independent constituent quarks the mass is of the order of 
N-Mgiuon- Hadrons with gluons can also be constructed {qqg, qqqg, ■ ■ ■)■ Of course, the 
gluons are highly virtual due to the strong interaction. It is not even known whether 
this picture is qualitatively correct. One can only refer to the quark sector and the suc- 
cess of the constituent quark model. Furthermore, a gluon mass modifies the transverse 
momentum distribution of gluon jets. With a theoretically computed gluon mass one 
has a natural energy cutoff. A gluon mass is also welcomed for instanton physics: In a 
massive gauge theory instantons of a size p > M^^^^ are exponentially suppressed. The 
infrared problem disappears. Under the assumption of the Instanton Liquid Model with 
independent instantons of radius p = 600MeV~^ and density n = {200MeVy one gets 
a gluon mass of Mgi^on = 480 MeV near the cutoff scale p~^ = 600MeV. Note that this 
possible solution of the infrared problem is different from the standard hope of finding a 



repulsive interaction between instantons [^, |2^, pO |. 



For a quantitative theory one has to generalize the calculation of unphysical gluon prop- 
agators to gauge invariant matrix elements, like glueball correlators. This is a problem in 
all works up to now (including this one). The methods to achieve a gluon mass are too 
complicated to be applied to physical quantities. For this reason it is difficult to make 
the above mentioned physical applications more quantitative. 

In the following Sections we will compute the averaged quantum gluon propagator for var- 
ious gauges in the multi-instanton background. Although explicit expressions are known 
for a long time [Q they have not been used in phenomenological applications since they 
are rather complex. Till now the instanton background has been treated at the classical 



level. Even the 1-instanton action ^Tj is of relatively low practical importance due to the 
infrared problem. On the other hand quantum corrections to gluonic n point functions 
can be large. 



7.2 Gluon Propagator 



* 



The first task to obtain a formal expression for the gluon propagator in a background 
field is to expand Cqcd[A + B] in the fiuctuations 5^ around our background A'^. The 
term quadratic in 5" is then by definition the inverse gluon propagator. For A'^ we will 
later use our instanton gas. We will use the QCD-Lagrangian 
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where we have rescaled the fields in such a way that the couphng-constant-dependence 
is in front of C We will entirely work in Euclidian space with metric Sf^u instead of (?^i, 
because instantons do not make sense in Minkowski space. At the very end we can simply 
rotate back to Minkowski space. With these conventions 

g'CQcoiA + B)= ^-G;M + B)G^,''{A + B) = 

0(BO) 0(b1) 0(b2) 



1 /^ — ''^ — ^ ' ^"^ — ^ 1 

4 



Gl^G^: + BlDfG^: +- i?;(-D-D;5,. - 2/,,,G^, + D^Df )i?^ + 



1 



^f,,,BlBlDfBt^-f^,,BlBl!,^,BiBl^d,,{..) , (7.1) 



where Z)^^ = d^bah + fach^p, is the covariant derivative with A inserted instead of A\ 
similarly G" = G" (A). As usual the terms which are total derivatives disappear after 
integrating L. We see that if A solves the equation of motion the linear term vanishes, 
as it should be. If we choose background gauge D^B^ = the last term of the 0{B'^)- 
contribution vanishes. Now we can read the inverse gluon propagator from the terms 
quadratic in B (from now on we will omit the bars over A, G and D because the unbared 
objects won't be needed further): 

iS-Xu = ^,i-D;^Df6,,-2UG;^) . (7.2) 

We will also omit the l/g"^ in front of the propagator which is a result of the rescaling of 
fields anyway. For further manipulations some abbreviations are useful: 



) 



{P^Y' = iDf , p, = td, , Pp = pp + Ap , 

PpX = [pp, X] + Xpp = t{dpX) + Xpp . (7.3) 

The last equation has only been quoted to show that p and P will be used in operator 
sense. F° are the generators in adjoint representation and fabc are the structure constants 
of the color gauge group SU{Nc). With these abbreviations we can now write 

S;, = P'S,, + 2zG^, = {p^+p.A + A-p + A^)6,, + 2tG,, = {S^' + V),, , 

-^J. = Wp' > Vp, = {A^+p-A + A-p)6p, + 2tGp, . (7.4) 

5*0 is the free gluon propagator with no background and V can be interpreted as an 
interaction potential caused by the background. The QCD-Lagrangian in the background 



7.1 can be written in another form more suitable for ordinary perturbation theory: 



g'CQCD = \G%iB)G^^^iB) + \B;v;tBl + UcfaeaB^BlBtA^^ (7.5) 
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The terms independent and linear in B have been omitted. Constant terms are irrelevant 
for the dynamics and terms linear in B are zero if A solves the QCD equations of motion. 
The new terms due to the background create two additional Feynman diagrams 




QjabcjaedA^d, 



up 



It should be noted that for small coupling constant g the second graph can be treated 
perturbatively but not the first. So our main concentration lies on the first term. 



7.3 Propagator in Statistical Background * 

The next step is to use some approximation scheme to calculate the propagator 

S = {S^^ + V)-^ = So{\ + T)-^ , T = VSo . (7.6) 

For large momentum p, Sq and therefore T are small and we can expand S in powers of 
T: 

S = So{l-T + T'^-T^ + ...) . (7.7) 

Note that S{x,y) = {x\S\y) is generally not invariant under translations and rotations 
because the background A^ and therefore T are not. Actually we are not interested in 
the propagator for a particular background configuration, but only in the average over all 
relevant configurations. Here we do not mean the functional integration over quantum 
fluctuations around the empty vacuum, but fields other than the perturbative ^4^ = 0- 
vacuum which minimize the action / C dx. 

5 = So(i-T + r2-r3 + ...) , (7.8) 

where the bar denotes averaging over relevant configurations. Details are specified below. 
If the background is statistically invariant under translations then S is translationally 
invariant and therefore diagonal in momentum space S{p,q) = {p\S\q) = S{p)^{p — q). If 
we would evaluate the terms in the series we would get an expansion of S in the form 

Sip) = l/p" + ci/p^ + C2// + . . . (7.9) 

and the pole remains at zero — and gets even worse with higher terms. What we want 
is S* in a form like S{p) = {p^ + M(p)^)~^ with M{p) bounded for large and small p and 
interpreted as momentum dependent gluon massQ. So let's invert 



S ^ =f + M{pf = (1 - T + T2 - T3 + . . .)-^Sq^ (7.10) 



^ The "+" in front of M will change to the more familiar "— " when we rotate back from Euclidian 
space to Minkowski space. 
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and expand it once again in T. Without averaging this would just be a geometrical series 
which, expanded, would give the original formula S^^ = Sq^ + V. With averaging the 
different terms are now in no relation and expanding and sorting with respect to powers 
of T yields 

S~^ = {l + T-{T^-T^) + {T^-TT^-T^T + T^))So^ + 0{T^) (7.11) 

= So' + V-(ySoV-VSoV) + ... 

= Sq^ + M"^, M"^ = Ml-Ml + ..., Ml = V, Ml = VSoV -VSoV. 

In the next section we introduce the concepts of instanton gas calculation to determine 
Ml, which is actually very simple. 



7.4 A Naive Estimate of the Gluon Mass * 



In the Instanton Liquid Model, in which A = J2i ^i is a sum of instantons in singular 
gauge with fixed radius p = 600MeV~^, the scatter amplitude 



Vf,^ = {A^+p-A + A- p)6^^ + 2iG^^ 
can be easily computed. A short calculation shows 



^M = (^/.)/ = , ^^. = (7.12) 

For Mf = F = ^2 we get 

CAi\ = 

For Ml = V = ^ this leads to 



(^ ) = 1^75 T'^P ^^b (7-13) 



Ml 



127r2/V 

-v^p^420MeV , (7.14) 



where we used the instanton density n = (200MeV)^. This value coincides very well with 



the results of 381 und 



For large momentum there is a term M2 which completely cancels Mi. A vanishing mass 
for large momentum is expected due to asymptotic freedom. To see this in detail one can 
count the number of p's occurring in M2 which will give us the dominant behaviour of 
M2 for large p. VSoV = Mf/p^ vanishes for large p but in principle V contains a p in the 
numerator {pA) and 5*0 = 1/p^ and VSqV could be finite for large p. To be more definite 
consider 

VSqV ={pA + Ap)Sq{pA + Ap) + other terms = AA^P^A^ + ... 

pi 
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where we have used [p, Aj] = id^A'^ = 0. 



{x\VSoV\y) = AA,{x){x\^\y)A,{y) + ... 

From S^u{x\Pf,Pu/p'^\y) = {x\y) = 6{x - y)we can conclude thai {x\pf,pi,/p'^\y) = \6^,y6{x - 
y) + traceless terms. 

{x\VSoV\y) = A^{x)A^{x)S{x - y) + . . . = {x\A^\y) + ... 

So M2 = A'^ + . . . = Mf + . . . contains a term which fully cancels M^ calculated above. 
Chapter |^ explains why such cancellations are expected in singular gauge. 

In the next section we present a systematic expansion of the propagator in the instanton 
density. 



7.5 Expansion in the Instanton Density 



* 



In this section we will make an expansion of the gluon propagator in the instanton density 
n which will be valid for all Euclidian momenta p especially for small p. Furthermore, the 
result will be Gauge invariant. 

The average of an arbitrary power of the 1-instanton field is proportional to the inverse 
volume ylj ~ p; for n > 1. So, for example, the expansion of the square of the multi- 
instanton configuration in the instanton density is 

N 

A2 = ^ A2 + iVA7 (A^ - 1)^7 = NAj +0{n^) (7.15) 

and more general 



0{n) 0(n) 0(n) 



A'' = NA^ + Oin^) , V" = NVp + 0{n^) for n > 1, (7.16) 

where V/ = V{Aj) defined in ( [7.4[ ) with A replaced by Aj and similar for T. Let us now 
sum up all terms in ( [7.11|) linear in n: 



5"' = {l + T-T^+T^-...)SQ^ + 0{n^) 

= {l + N{T;-Tf + Tf-...))S^' + 0{n^) 

= {l + NT:ff)So' + 0{n^) = So' + NV:ff + 0{n') (7.17) 



Teff = Ti-Tf + Tf-... = Tj-Ti{Tj-Tf + ...)= Ti-TiTeff 

Veff = Vj-VjSoV,ff =^ V,ff = S,\So-Sj)S,' with (7.18) 

5*7^ = Sq^ + Vi is the propagator in the 1-instanton background. (7.19) 
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More generally we can make a cluster expansion of an arbitrary function of A in the 
following form 

N 



f{A, + ... + A^) = f{A) = E(f ) Y: i-y-'QfiA, + ... + A,)= (7.20) 

1=0 k=0..l 



= /(O) + NfiA,) - /(O) +- NiN - l)/(Ai + A2) - 2/(^2) + /(O) +Oin') 

"— v—^ ^^ V ' 2 V ' 

0(1) 0{n) Oin'^) 

where the first line is an identity even without averaging. Inserting a Taylor expansion 
for / in the second line and using the indistinguishability of different instantons one can 
see that all monomials in the k-th term contain more or equal than k different instantons. 
So the average will factorize in k factors each proportional to n and therefore the k-th 
term is indeed proportional to n^. It is easy to generalize ( [7.20| ) for two or more species 
of fields. If A is a field of Nj instantons and Nj anti-instantons we get to first order in n 



f{A) = /(O) + NjfiAj) - /(O) + Nj-f{Ar) - /(O) + 0{n'). (7.21) 

If we insert the propagator S in (|7.21| ) we get 



S = S{A) = S{0) + NS{Ai) - S{0) + 0{n') = So + NSj - So+ Gin') (7.22) 

which is after inversion up to Oln^) just ( |7.17| ). 



7.6 QCD Propagators * 

In Section 5 we have seen that it is enough to know the 1-instanton propagator to calculate 
S to first order in n. Luckily this propagator is known, unfortunately it is a bit lengthy 
and suffers from a divergence. 

The gluon propagator Sf^^ with spin S* = 1 in adjoint colorQ representation (C = 1) can 
be constructed out of the ghost propagator Aj* {S = 0,C = 1) which is explicitly known 
in the 1-instanton background. The general formulas how to construct a propagator of 
given spin S out of the corresponding scalar propagator with same color C in a selfdual 
background are shown below. They are derived and more thoroughly discussed in |TB|. 



Notations: 



A^ = T'^Al , [T^T^] = ^e„,eT'= , T« = generator of SU{2), 

in scalar (C = 0) 

rpa ^ ^ ^a/2 in fundamental {C = i) 

ie.a- in adjoint {C = 1) 



> representation (7.23) 



-Pa* = Pm + ^M ' Pm = ^^M ' ^fJ-i^ = T^^i^upcrGpa , {lpi,lu} = 2(5, 



2 



-'fiu 



^ Often called isospin pi 
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Spin propagator A : 

A"^ = p2 =^ A = p-2 (7.24) 

Spin I propagator S : (not used but only stated for completeness) 

S-^ = If=^^P^ =^ 5 = ^Ai^ + Afi^ for G,, = G,, (7.25) 

Spin 1 Propagator S: 

S-^ = P'A^, + 2tG,, - (1 - hp^P, =^ 

S,, = q,,p„PpK^P,-{l-OPA^Pu for 6-^. = ^^., (7.26) 

There are some comments in order. The above formulas are vahd for an arbitrary color 
representation, we will need them only for C = 1. For S ^ 1 there are zeromodes and the 
propagator is only the inverse of the kernel in a subspace orthogonal to the zeromodes. 
The implications and problems will be discussed in Section 8 when they show up explicitly. 
The Spin 1 kernel is the quadratic term of the QCD-Langrangian ( |7.1| ) with gauge fixing 
term ^(L'^'^S^)^ in slight generalization to the ^ = 1 case considered previously. 



With 



n(x) = l + 4 , F{x,y) = l + p^^^^^ , (7.27) 

X X y 



r„ = [T,i 



p \' 1 "I 1 ' p 



rl = (r, -i) , t^tI = 5^^ + irfapuTa 



we can write the ghost propagator for 1 instanton in the form |]T8 
kr TaF{x,y)nF{y,x) 



Af{x,y) 



47r2(x-y)2n(s)n(y) 

Sab p'^Sab '^p'^eabcVcpuXpyu 



+ 



An'^[x — yY 47r2(x^ + p^)(?/2 + p^) Ati'^[x — yY^x"^ + p'^){y'^ + p'^ 

2p'^{{{xyY - x'^y'^)5ab + tabcVcpuX^yy + f]apuVbpaXpyuXpya 



47r2(x — yY{x'^ + p'^)x'^{y'^ + p^)y^ 



(7.28) 



For simplicity we have placed the instanton at the origin z/ = in standard orientation 
Qab _ ^ab j^ jg possiblc to Write down the gluon propagator S'}^ explicitly in which we 
are finally interested in using ( [7.26|) but the expression will be rather lengthy. While /S.'f 
can be averaged and represented in momentum space exactly with the help of modified 
Bessel functions this seems not to be possible for Sfl^^^ and we have to expand S for large 
and/or small momenta (or work much harder to solve the complicated integrals in terms 
of special functions). So we will never use the full expression for S. 
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7.7 Propagators for Small Momentum * 

The calculation simplifies significantly if we expand A{p) or S{p) for small momentum. 
Because p always occurs in the dimensionless quantity (pp) the lowest order in p can be 
found by keeping only terms of lowest order in p or differently stated: Small p corresponds 
to large x and for a; S> p p is negligible. To warm up let us start with A from ( |7.28| ): 

Af = A;^'' - p^W'' + Oip^) (7.29) 



,abf^ „\ ^ah 



Ati'^^x — hY 



rabf \ o-b ^^abc^c^u-^ ^Vp 



W^"(^,^) = 7r£b^ + 



Ao(p) = 1/p^ is the free (A" = 0) ghost propagator. After reintroducing the instanton 
position z we can now average Aq — A/. The e-term in W will be killed by averaging over 
the instanton orientation: 

(x|A?rA?b) = 4.^W^l/j^,^^^-i^^,i^-^ . (7^30) 

The last integral is infrared divergent but noticing that l/47r^(a; — y)"^ is the Fourier 
transformation of 1/p^ we can rewrite the above expression in the form 

V^{x\Ao-Aj\y) = 47rV f d''z{x\\\z){z\\\y) = 47rV(x|^|y) 

J p'^ p'^ p* 

The divergence is now hidden in the fact that the coordinate representation of p~^ does not 
exist. In fact we are not interested in the coordinate representation but in the momentum 
representation and so the divergence is spurious. Inserting A in ( [7.18|) we get 



Vl)f = A,' A^^AT A,' = P'^^y = ^^^^y = pYW (7.31) 

and lead in the SU{2)c case to a ghost mass of 



2 „2„ 



Mt,,ost{P = 0) = NV^ff = A^'p 



n 



which has to be inserted into the ghost propagator A^^ = 5ah{p^ + ^Ihostip)) ^ ■ General- 
ization to the case of N^. colors gives 

Mghost{p = 0) = 27rpv/2^ = 8.9pv/^ = 340 MeV (7.32) 

where ur := n/Nc is of the order oiO{N^). The most important insight is that M ^ 0. A 
scalar particle in adjunct color representation achieves a dynamical mass in the instanton 
background. 
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Let us now calculate the gluon mass at zero niomentuin along the same lines. To do this 
we must insert (|7.29| ) into (|7.26|) . Expanding also P^ up to order p^ 



P,=P, + p'A, + 0{p') , A, = ^^"^'°f ^- 

we get 

Pf^AJP, = vAoVu - p'v^.{^oW + iyAo)p. + p'(p^Aoi. + A^A^p,) + 0{p'') 

The free gluon propagator in i?^-gauge is well known or can be obtained from (|7.26|) 
setting G^u = 0: 



S' - S' = p'[q,up.Pp{AoW + WAo)p^ - (1 - OP.i^oW + WAo)p, 



S% - Sf,, = 2p'W{p%, - (1 - Op,p.) = ^pVwsl 
where we have used in the last line that p commutes with averaged quantities like W. 



p'So V-(l-e)^ 

We got the interesting result that for momentum transfer the gluon mass is a factor of 
a/2 larger than the ghost mass. Maybe the relation Mgiuon ~ ^Mghost is valid even for 
larger p. The gluon mass for momentum transfer zero is 

Mgi^onip = 0) = 47rpy^ = 12.6pv^ = 480 MeV (7.33) 



7.8 Zeromodes * 

In our whole calculation we have ignored the zeromodes. All gluon field fluctuations 
must be orthogonal to these which is achieved by using a gluon propagator orthogonal to 
the zeromodes or otherwise stated by subtracting out from ( [7.26| ) the projections on the 
zeromode subspace. This procedure also deletes a divergence coming from the second term 
in ( |7. 28| ) squared. There are some further problems to be solved because the scalar product 



of the zeromodes with the propagator does not exist |T8[. But all this concerns a term 
which is proportional to p^ which was irrelevant in our zero momentum approximation. 
Also the principle mechanism of mass generation does not depend on the zeromodes 
which can be seen from the ghost propagator because in the spin case there are no 
zeromodes. This should be contrasted with the fermionic case where the zeromode are the 
main ingredient for mass generation and the so called zeromode approximation simplifies 
calculations enormously. 
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7.9 Conclusions and further developments 

We have computed the gluon propagator in the Instanton Liquid Model for small mo- 
mentum p in leading order in the instanton density. A mass independent of Nc has been 
extracted. Furthermore, the mass is gauge independent, at least for o^ = within the 
class of R^ gauges. As expected, for ^ — i> 0, the gluon propagator is transversal. 

A massless pole remains for ^ 7^ 1 due to the (1 — p)p^Pu/p'^ term. We expect that 
the massless state, which is associated with this pole, decouples from physical matrix 
elements, as in conventional perturbation theory. 

The result is in agreement with the values calculated by Cornwall (500 ± 200Me\^) 
or extracted from pp scattering (370Mel^) 



The next step may be the calculation of some gauge invariant correlation function like 
a glueball correlator or topological susceptibility. As usual this involves products and 
integrals over propagators but with the difference that we have to use ( |7.26| ) instead of 
the simple free propagator. Without further simplification this may cause a headache. 



Chapter 8 

Gauge Invariant Quark Propagator 



A variety of predictions concerning chiral symmetry breaking and concerning the lightest 
hadrons in various channels can be made within the instanton liquid model. Although 
there are various attempts to derive this model from first principles, it is still an open 
question, whether instantons melt or not. Thus the infrared problem remains unsolved. 
In this chapter I attempt to make a few predictions which do not rely on the Instanton 
Liquid Model, but exploit the theoretical one- loop density D{p). 



The quark condensate will be calculated in Section ^]5| and p^ . I get a finite result 
by choosing an appropriate gauge and performing a self energy resummation without 
introducing an infrared cutoff in the instanton density. 

Because the finiteness essentially depends on the choice of gauge, I give a more general 
discussion in Section JO], ^]2] and |8]^ of how to choose a gauge when calculating gauge 
dependent quantities. The quark propagator in the background of one instanton in the 
well known regular and singular gauge is compared to a gauge invariant propagator, for 
which explicit expressions are calculated in Section |8.1| . 



8.1 Generalities On the Choice of Gauge 

Gauge symmetry is a rather large symmetry, an infinite product of SU{Nc) in the case of 
QCD. A physicist is always happy of having symmetries because they can be exploited to 
make predictions even without solving the theory. Gauge symmetry is necessary to get a 
physical vector particle spectrum. As long as one does not make an approximation which 
manifestly breaks gauge symmetry one can choose a comfortable gauge for calculations 
because the result is GI. But it is very difficult not to break GI, especially in a non- 
abelian gauge theory. It is not easy to find a GI regularization and furthermore, the gluon 
propagator, the primary object in perturbation theory, is not GI. Of course it is meanwhile 
well known how to perform GI calculations in every order perturbation theory using FP- 
ghosts and dimensional regularization. Every new approach beyond perturbation theory 
is again confronted with the problem of GI. In lattice theory the Wilson action had to 
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be invented. In Schwinger-Dyson and Bethe-Salpeter type selfconsistency equations GI 
is still an open problem. In instanton physics when going beyond the one instanton 
approximation the choice of gauge is also important. This will be discussed in the next 
paragraph. There is a related problem when considering non-GI objects from the very 
beginning like the gluon or quark propagator. Strictly speaking they are only defined when 
relying to a certain gauge. In principle one should not give them any physical meaning at 
all. Often one is tempted to do so and therefore it is necessary to give some motivation 
of choosing this or that gauge. 

8.2 A Natural Gauge 

The gauge field A" describes the connection between neighboring vector bundles over the 
spacetime manifold IR^. Thus a choice of gauge is like the choice of a coordinate system 
in general relativity with connection F'^ . When choosing a crooked coordinate system, 
although being in a smooth universe, there will appear fictitious accelerations towering 
above the real physical accelerations 

•'^phys ~ -^fict "T ^ vp-^ -^ 

When making general covariant calculations these fictitious accelerations and the F con- 
tribution will cancel out thus leading to the correct small result. But the slightest un- 
systematic approximation will produce gross errors. The natural solution of this problem 
is to use a coordinate system as smooth as possible to avoid fictitious accelerations, e.g. 
to choose F'^ as small as possible. To make this statement more quantitative we may 
try to minimize (x^^y^ — x^fictY simultaneously for all curves. This is done by choosing a 
coordinate system which minimize^ 

llriP •= [ V "T^ d^r 

This obviously measures the crookedness of the coordinate system. 

Let us now transfer this to QCD. The analog norm for the gauge potential is 

PIP := Jtr.A.A^'d'x 

A stationary point is found by variating \\A\\ w.r.t. gauge transformations 

6Af, = i[A^,Q]+df,Q , 6\\A\\^ = 2t f tr{d^A>')Qd^x = Vfi ^^ d^A" = 

Therefore in Lorentz gauge, A"^ contains as few pure gauge as possible, if the stationary 
point is a minimum. An expansion in A is thus most rapidly convergent in Lorentz gauge. 
In applications where A is not needed in total e.g. when only a certain momentum region 
is probed, different norms and different gauges may be optimal in the sense discussed 
above. Especially one should include derivatives of A into the norm in order to guarantee 
a smooth A which is important for high energies. 

^In Euclidian space this is a positive definite norm 
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8.3 On the Gauge in Instanton Physics 

When calculating GI quantities in the background of one instanton in a GI way the choice 
of gauge is only a matter of convenience. But one can see that there are large cancellations 
between different terms in regular gauge at large distances due to their slow decay and in 
singular gauge at small distances due to the topological singularity at the instanton center. 
For non-GI invariant quantities like the gluon or quark propagator, or when making some 
unsystematic approximation, the lesson is to use singular/regular gauge when dealing 
with low/high energies to avoid these cancellations. This is consistent with the discussion 
given above. Singular as well as regular gauge fulfill the Lorentz condition. ||y4sj„g|| is 
finite and a minimum. Agi^g is therefore a good choice for low energies. For high energies 
it is important to have a smooth A which is obviously only satisfied by the regular gauge. 

To linearly superpose instantons they have to decay rapidly enough. Therefore one has 
to use singular gauge. This argument can in principle be circumvented by superposing 
two fields An and Ajj- the former/latter being an exact multi-instanton/anti-instanton 
configuration in regular gauge. Despite this, for low energies singular gauge is in any 
case a good choice and for high energies a one instanton approximation is already a good 
approximation. 

8.4 The Quark Propagator in Axial Gauge * 

A specific example to test the gauge dependence is the quark propagator. The contribu- 
tion of one instanton of radius p to M{p) = ip'^Sj{p) which usually is interpreted as a 
constituent quark mass is shown in figure [B.2| in regular, singular and axial gauge. The 



regular graph is larger than the singular at low momentum and the singular graph shows 
the slow decay (only polynomial in 1/p) for large momenta. The analytical expressions 
are well known and are listed in Appendix |A.2| together with expressions in axial gauge 
which will be derived and discussed below. 

A correlator containing color- non- singlet operators can be made GI by connecting distant 
points with a special path-ordered exponential containing the gauge field. The exponen- 
tial ensures the parallel transport of color from one point to the other. The GI quark 
propagator may symbolically written as 

SaA^.y) = (0|^(a;)Pexp (^t J^dz-Aiz)^ ^{y)\0) (8.1) 

P denotes path-ordering. We have already defined Sax to be its color singlet part because 
only the singlet part is GI. Sax will be called the axial propagator because in axial gauge 
with n^ = x^j_ — y^ the exponential vanishes. In the one instanton background in zeromode 
approximation we get 

Sax{x,y) = j^tTc 



Pexp ( i / dz-A{z) ) ip{x)ip{y) 
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where A is now the instanton field and ip is the zeromode in any gauge. In a coordinate 
system where the instanton sits at the origin and x — y is in time direction^ (x = y = z) 
the path ordered exponential reduces to an ordinary exponential. Alternatively we could 
have tried to find a gauge transformation which transforms the regular gauge in axial 
gauge. In both cases we get: 

Sax{x,y) = -^tr^ i)ax{x)i)ax{y) , i)ax{x) = R{x)i)reg{x), (8.2) 



N, 



c 



R{x) = e *°^^^T^ = cosa(x) ± i- — - sina(x) =: ±ir^-x(a;) (8.3) 



/ \ "^ / ^ \ 

aix) = . ,. ,r arctan . ,. ,. : (8.4) 

Vx^ + p^ Vx^ + p^ 

a{x) may also be written in a covariant form 

( p\x-yf \~'/' 
a[x) = ± 1 + -^-^ — — arctan 



(a;2 - {xy)Y 



x^y"^ — {xyY j \ x^y^ — {xyY + p'^{x — yY 

but now a{x) depends also on y and the expression for the propagator no longer factorizes. 
The reason for this is that the axial gauge is not covariant, but the definition of the 
propagator is. Inserting ( ^.3|) and ( ^.8|) into (|8.2| ) we get 



Sax{x,y) = — Uxy) - -Xf.y^a'"'] ^(firegix)ipreg{y) (8.5) 

Inserting (p.3|), ( |8^ ) and ( p.8[ ) into ( |8.5| ) the space-time averaged propagator can be 
expressed as an integral over elementary functions 

Fi/ \ f^ f°°, A 2 \^ ( t+\x-y\ t^ 

b{x — y) = dr at Anr cos — arctan arctan — 

JO J—oo -^ V -^ -^y 

1 o^ 

R = r + p 



2Nc 7r2(/?2 + (t + |a; _ ^|)2)3/2(^2 + ^2)3/2 

The difference between the propagator in regular and axial gauge is the insertion of the 
cos[. . .] factor. Therefore the axial propagator is everywhere smaller than the regular 
propagator, except at x = y where they coincide because the path-ordered exponential 
is one. At large distances it is smaller by a factor tt/4. Instead of performing the inte- 
gration in coordinate space, let us go directly to the more interesting momentum space 
representation: 

5,(P) = ^^^.(p)^:^Kp) , Vi^AP) = J i'{^)¥>reg{x)e^'''dx (8.6) 

Although (y9^^ does not transform like a vector, we can choose a convenient direction of p 
because ifip"^ is a Lorentz scalar. For pure spacelike p the spatial components of ip vanish 

^ Although working in Euclidian space we will adopt the Minkowskian language (xo,x) — (time, 
space) . 
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because the integrand is anti-symmetric w.r.t time reflection. Only the time component 
is nontrivial 



With the following hints 



d\ / dt cos 



r t ' 

— arctan — 
R R. 



P 



7r(i?2+t2)3/2 



o«pr 



/oo 
{R 
-oo 



cos(7 arctan x) 
ity{R + it)-^dt 



Re 



1 + ix\ 
1 — ixj 



7/2 



2n{2RY 
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the reader should be able to perform the t and the angular integration ^^2^ 

/vrr \ 



•7) 



y^axip) 



ppJo 



oo 
COS 



] sm(pr)r dr 
\2RJ 



I was not able to perform this last integral analytically, but for small momenta it is easy to 
see that v^a^-d^) behaves like ir'^p/p. For large p it decays like ~ e~^'' with a non-polynomial 
coefficient because of an essential singularity at r = ±ip. Comparing ipregi fsing and ipax 



plotted in figure p.2| we see that the axial (p lies somewhat in between the regular and the 
singular. So one may conclude that axial gauge is a good compromise for all momenta. 

The calculation of the GI propagator seems to make the discussion of its gauge dependent 
partners obsolete. I will now argue that this is not the case. The reason is that there 
are a huge number of GI definitions of a quark propagator and ( p.l|) is only one possible 
choice. One obvious generalization is to choose a more complicated path from x to y than 
a straight line. The next thing one could do is not to restrict oneself to a specific path, but 
to take into account all paths one is interested in and average the results with arbitrary 
weights. Another possibility is to let the path depend on the gauge field itself, as long as 
this choice is made in a GI way. Finally one can combine both generalizations. I am sure 
that it is possible to produce any result for the propagator with a suitable generalized 
definition. The advantage of the standard axial propagator is, that the definition is simple 
and that the non local operator has a physical interpretation. It creates a quark-antiquark 
pair connected by a thin gluon flux tube. This might be a good choice for a non-local 
meson creation operator. But it is also plausible that one of the generalizations given 
above is even better. The only thing I want to point out is, that the GI definition for 
the propagator given above is nothing more than to work in axial gauge. One still has to 
choose the right gauge using more sophisticated arguments. 



8.5 Effective Quark Mass * 

The only reason for performing the N^ ^ oo limit is to make the A''^ dependence of the 
resulting formulas simple. The accuracy has been checked to be within the standard 
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10% for Nc = 3 usually achieved by 1/Nc expansion. Here the accuracy can simply be 
understood. The actual expansion parameter is not 1/Nc = 1/3 itself but 1/b ~ 1/11. 
The following asymptotic formulas will be used 

iVji/^^ ^ iV,/e , b = —N, , C^(' ^ 2.22r^/" 

p'D{p)r.{2.22{So/bf/''pAf , So/b='^{g',Nr' (8.9) 

Every equality in the large Nc limit will be marked with a dot. Notice that in this limit 
instantons of size p < ■^{So/b)~^^^^A~^ are completely suppressed. Above this threshold 
the instanton density gets infinite. So/b is independent of Nc because the coupling is 
go ~ l/^/iVc- A is the QCD scale. 

In the presence of one light quark flavor the instanton density D{p) has to be multiplied 
with the functional determinant of the Dirac operator 

Det{ip+ im) ^ l.SAmp 

which is proportional to m because of a zeromode of p. The quark propagator in the 
background of one instanton is dominated by this zeromode 

im 
Averaging this expression over all collective coordinates 7/ one gets 

M{p) ■= tp^Siip) = -— dpp^pD{p)ip\p) 
2i\c Jo 

Summing the contribution to the propagator of 0, 1,2,3,... instantons, which is the analog 
of a selfenergy resummation in perturbation theory, 

S(p) = 1 + l^^l + 1 M(p^ 1 M(p^ 1 ^ ^ 1 

^ j^ i j/ 1/ i 1/ i 1/ " ' ll + iM{p) 

justifies to call M{p) a dynamical quark mass. Expressions of (p in various gauges are 
given in appendix |A.2| . The graphs of (p{p) in singular and regular gauge cross over at 

'Psingip) = '^reg{p) ^^ 2pp ^ 2.5 

Therefore one should use regular gauge for large p and singular gauge for small p. This 
choice of gauge also makes the integral convergent for large p. At this stage we have no 
infrared problem. Using regular gauge in the whole integration interval we get 

MreM = Bpi^f , B = 1.34 ■ 167r\CNy''^){So/b)''/''h/Nc 
zp 

CO 

dzz^~^e-' = {b-2)\ , z = 2pp 
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The integral is sharply dominated by 

z = 6±6V2 »2.5 

therefore the result is independent of the choice of gauge for z < 2.5 justifying our use of 
regular gauge over the whole integration interval. Axial gauge would lead to nearly the 
same result as can be seen from figure |B.2| . Using singular gauge for large p would produce 
a divergent integral dominated by arbitrary large instantons inconsistent with the choice 
of gauge discussed above. The infrared "problem" shows up in the rapid raise of M{p) 
for low p, which effectively suppresses the propagation of quarks with low virtuality p"^. 
Consider some process involving quarks at distances x = 1/pc- The effective quarkmass 
M{l/x) is dominated by instantons of much larger size 

p = Pe(l±rl/2) »X , Pc = — • 

In other words, given an instanton of radius p influences the physics at a much smaller 
scale X = ^p << p. Therefore the interior of the instanton is probed and one should avoid 
the singularity at its center by using regular gauge. 

8.6 The Quark Condensate 

Let us now calculate a real physical gauge invariant observable, the quark condensate 

(m := limt.co(S(.) - So(.)) = -4«o/ ^j^fLj^ 
Inserting M(p) and performing the angular integration we get 



J, = / ^±-^dz = -.— 4n^- ^ 4 , p = B'/'-z 



l + 22b 26sin(^7r) 2h ' ^ 2' 

The integral is finite and sharply dominated hy z = 1 ± h^^ . Without resummation of 
the selfenergies the integral J^, and thus condensate would have turned out to be infinite. 
The condensate is dominated by quark wavefunctions with momenta 

P = Pc{l±b-^) , Vc = Ph^ , /3 := >/^ = ^(V&f/" 

2 be 

and depends on A and qq 

|(7/;7/^)|i/3 = 0.139/36A . 

Expressing p^ and pc in terms of | {^ip) \ by eliminating (3 we get our main result 

Pc = 3.59\{ijij)\^/^ 

Pc = ^l(V5^)r/^ (8.10) 

2.22 Ar,A = {g',N,f/''\{ijij)\'/' 
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A weak point is the experimental extraction of go- It should be extracted from a reliable 
treelevel process at low energies presumably of the order of pc- In QCD improved Bag- 
Models the main nonperturbative effect is modeled by the bag and the hyperfinesplitting 
is caused by a one gluon exchange, go extracted from A — A^ splitting is |10| 

go"' ^ 2.6 

Let me also give a theoretical guess of go. The change to a two loop expression for the 
instanton density So/i '^ S1/2 can be effectively performed by only replacing So in the 
following way 

Because a is very small (In ^)" is approximately one in a large range of values for pA. 
and for 



the two 2 loop density coincides with the one loop density. Because we do not believe 
that the 2 loop density is an improvement, one should not take 5'o"'^'^^ too seriously. At 
least it is not in contradiction with go""^. 

The condensate is well known to be \{ijjip)\^^^ = 240MeV. Setting A^c = 3 and taking 
(70 = 2.6 for granted we get 

Pc±Ap = (860±80)MeV 

Pc±Ap = (160 ± 50)-^MeV (8.11) 

Apv ^ 190MeV 

The most interesting thing is, that the condensate is sharply dominated by quark field 
wave functions of rather large momentum pc- On the other hand the dominating instan- 
tons have a very large radius pc, 4 times larger than usually assumed in instanton liquid 
models. Nevertheless the predicted value of Apy, which of course must be assigned a large 
error because of the rough estimate of go, is in agreement with experiment. 



8.7 Summary 

Whenever one is calculating gauge dependent objects or when making gauge breaking 
approximations, one is confronted with the problem of choosing a "good" gauge. Spe- 
cializing the general discussion of Section ^.2| to the case of instantons, we came to the 
conclusion that the regular gauge is appropriate for small distances and the singular gauge 
for processes involving large distances. The GI propagator was defined, calculated and 
compared to the propagator in singular and regular gauge (figure |B.2| ). The conclusion 
was, that the GI propagator is not a-priori a good choice, but lies somewhat in between 
regular and singular gauge. 



Using an appropriate gauge along the lines discussed in Section |]2| we were able to derive 



a finite quark condensate without taking an infrared cut-off for the instanton radius nor 
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relying on some instanton model. The linear relation between \{^jjip)\^^^ and the QCD 
scale A is in agreement with experiment. The condensate if formed by quark fields of 
high momenta pc = 860MeV mainly lying within the sharp region Ap = SOMeV. The 
dominating instantons are very large {pc = 160MeV). 



Chapter 9 



Conclusions 



9.1 Summary 

Some known and various new results have been obtained in this work. In most cases, the 
computation was based on the Instanton Liquid ModeL 

The following new insights were gained. 

• For the quark propagator dynamical quark loops can be absorbed in a renormalized 
instanton density, which can be identified with the gluon condensate (Chapter]^). 

• Concise Bethe-Salpeter equations have been obtained and were solved for the quark 
4-point functions. In the fiavor singlet channel, a chain of quark loops contribute, 
which is responsible for the missing f/(l) Goldstone bosons (Chapter U). 

• From the meson correlators the masses and couplings of the cr, p, cj, ai and /i 
mesons have been obtained by a spectral fit (Table |5.3| , Figures |B.3| - p.8| ). 



The Tj' mass has been predicted successfully (Equation |6.16| ). 

The axial proton form factors of the singlet current j^s, the gluon current i^^, 
and the anomaly a have been computed and their gauge (in) dependence has been 
discussed (Table |6^ ). It has been shown that ^(0) = —1 is independent of the 
number of fiavors. 

For small momentum a ghost and a gauge independent gluon mass have been cal- 
culated (Equations [7.32| and |7.33| ). 



General rules regarding the choice of a gauge, especially when to chose the singular 
and when the regular gauge have been established (Sections |8]^, |8]^ and |0| ). 



A gauge invariant quark propagator in the 1-instanton background has been derived 



(Equations ^j^ and p.8|) . 



The quark condensate has been identified with the QCD scale A, where neither an 
infrared cutoff, nor a specific instanton model was necessary (Equation |8.10| ). 
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9.2 Outlook 

The Instanton Liquid Model has proven to be successful in various sectors of QCD as the 
results of this and other works show. Although this model cannot strictly be derived from 
first principles, its success proves that there is some truth to it after all. The limits of 
this model also show up clearly: Confinement cannot be explained and the axial singlet 
channel still causes problems, although instantons explicitly break the U{1)a symmetry. 

Future work should take into account non-zeromodes, which become important in the case 
of strange quarks and probably essential in respecting the axial Ward identities. Baryon 
correlators to determine Baryon masses could be calculated. Both extensions have already 
been realized numerically p8| . 



Although a direct calculation of the axial singlet correlators was not successful, refined 
arguments allowed determining the t]' mass. A similar consideration may also solve the 
proton spin problem. 

Gauge invariant glueball corrections should be calculated. 

The most urgent theoretical problem is still a correctness proof of the Instanton Liquid 
Model. This would clarify the range of validity and increase the general acceptance of 
this model. 

QCD is, indeed, a tough theory. Experimental physicists do a good job in measuring the 
hadron parameters, whereby the precision of these experiments increases steadily. This 
progress should challenge the theorist to keep up with them. 
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Appendix A 



A.l Notations 

Nearly all occurring factors 2tc can be absorbed in the following definitions: 



{2nY ' 

At many places vectors and operators in Dirac notation in IR^ are used: 
{x\p) = e'P'-" , {x\y) = 6\x - y) , {p\q) = t{p - q) 

{x\ij)=i){x) , {p\-^)=-^{p) , {p\S\q) = S{p,q) 
d^p\p){p\ = l. , d^\x){x\ = l , 

[p,X]=tid,X) 

This notation should not be confused with the vacuum state |0) and the proton state \ps) 
of the QCD Hilbert space and the averaging (■ ■ ■)/ over instantons. 



iVc = 


= Number of colors 


Nf - 


= Number of quark flavors 


m = 


= Current masses 


M = 


= Dynamic masses 


tlD = 


= Trace in Dirac space 


trc = 


= Trace in color space 


Tr = 


= Functional trace 


Det = 


= Functional determinant 


A72 = 


= Generators of SU{N^), a = 1 . . . N^ - I 


r72 = 


= A72 = Generators of SU{2)c, a = l... 
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A. 2 Instantons in Singular, Regular and Axial Gauge 

Instantons are solutions of the classical Yang-Mills equations of motion. The instanton 
at the origin in standard orientation in singular, regular and axial gauge is given by: 

2 

A;^^(x) = R{x)A'^^^{x)R\x)+iR{x)df,R\x) (A.l) 

The upper/lower sign denotes an instanton/antiinstanton {Q = ±1). 

-I- „ / N „/ N / cosa(x) \ 

R{x) = ±iT^x^{x) , x^x =( . \\\ 

a(a;) = arctan ^=^=^ 

The covariant derivative p has a zeromode 

where the zeromode has the following form: 

iJregix) = ^Plif{x)x , V{x) = ^^^^ /^2)3/2 (^-2) 

^ax{x) = V2(p{x)R{x)x 
X is a color & Dirac spinor, given by 



16 



-^ -^ ~ 1 c ^M^i' 9 '''^ ''~i^ 



For light quarks the propagator is dominated by the zeromode. When averaged over the 
instanton orientation, position and charge the propagator is diagonal in momentum space 
and given by 

where 

i'singi.P) = V2v^sms(p)-r-|X , (fisingip) = T^P^ -7-[h{z)Ki{z) - Io{z)Ko{z)]^=pp/2 
tpregip) = V2(pip)x , ^regip) = e'^'' (A.3) 

8 f°° f nr \ 
i^ax{p) = FT{iJax}{p) , ^ax(p) = — y^cosl — ==^ lsin(pr)rrfr 



A.3. AVERAGING OVER THE INSTANTON PARAMETER -^ 
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l^{p) 


singular 


regular 


axial 


pp<^ 1 


27r 


4:71 


vr^ 


pp> 1 


1277 


Aire-PP 


~ e~Pf 



Table A.l: Asymptotic behaviour of -^p{p) 



Iv{z) and Ky{z) are modified Bessel functions. The asymptotics are given in the following 
table 

The constituent mass of a quark in the Instanton Liquid Model is 

Mp(p)~pV(p) ■ 

Mp{p) is plotted in figure |B.2| in all three gauges. 



A. 3 Averaging over the Instanton Parameter 7/ 

The instanton in general orientation and location and the corresponding zeromode have 
the following form 

Ai^{x) = OfAl{x-zi) 

iPiix) = Uj^{x-zj) , Xi-=UiX , ^tr(f/A<^f/tA^) = O"^^ 
7/ = {zi,Oi, pi,Qi) = (location, orientation, radius, topological charge) 

A^ is and instanton in standard orientation with radius p = pi and charge Qi = ±1 



and ip is the corresponding zeromode, given in Appendix |A.2| . Uj G SU{Nc) and 
Oj G Ad[SU{Nc)] are orientation matrices in fundamental and adjunct representation, 
respectively. 



t various points I have to average over the collective coordinates 77: 

(. . .)/ = fd^j D{pj) ... = I Yl fd%dOjdpj D{pi) . . . 

The following assumption on D defines the Instanton Liquid Model: 

B{pi) = n5{pi - p) , n = (200 MeV)^ , p = 600 MeV"^ 
The most important formulas for the Haar measure / dOi and / dUi are: 

I dO 1 = 1 , I dO O^^O"^ = j^2^_J ""5'"^ 

J du 1 = 1 , J du uiul^ = j^sisi 



(A.4) 



(A.5) 



(A.6) 
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Integrals over an odd number of matrices are zero. The following formulas are also useful 

Nc{xfxf)i = trcxfxf = ^(^^) (A-r) 

where the averaging over Qj has not been performed yet. 

XX = trcDXX = 1 (A.8) 

Collecting space, color and Dirac terms yields 

4ip)Mp) = '^f\p) , J(^4{p)Mp) = '^ (A.9) 

A. 4 Numerical Evaluation of Integrals 

The integral expressions for the meson correlators have been evaluated numerically. Two 
types of operations have to be performed: 

1. Convolution of Lorentz covariant functions (Fo/5,rr) 

2. Fourier transformation (FT) of the correlators to coordinate space 

Let us first consider the FT generalized to d dimensions: 

L...,M = F,{f,,,„,J{x) = /^e-^V,....M„(p) (A.IO) 



For a scalar spherically symmetric function / = /(|p|) the FT reduces to a one dimensional 
integral 

Fa{f{\p\)}{x)= i-y-] f{m)J,/2-i{m\x\)\x\dm (A.ll) 

where Ji, are Bessel functions. 

If X is not too large and / decays rapidly, the integration can be performed with Gaussian 
(or other) integration methods. If the decay is too slow one has to subtract the asymptotic 
part from / thus improving convergence. The FT of the asymptotic part can be performed 
analytically and has to be added to the numerical FT of the reduced function. 



The FT of a general Lorentz covariant function can also be reduced to ( [A.ll|) with (for- 
mally) an increased dimension d: 

F,{pJ{\p\)}ix) = 2mx,F,^2{mx) (A.12) 



1 



J^iiJy 



Fd{p,pJ{\p\mx) = j^[{S,.~^)F,{p'f{\p\)}{x)- 

(V - rf^)(4vrF,+2{/}(x) - 47rVF,+4{/}(x)) 
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A. 5 Numerical Evaluation of the Convolution 

The convolution integral is defined as 

f*g{p) = jc^qf{q)-g{p-q) (A.13) 

This type of integral can be reduced to the FT discussed above: 

f*g{p) = F,^{F,{f}.F,{g}} (A.14) 

This is a quick and easy method for evaluating convolution integrals. For the disconnected 
part of the correlators it has the advantage that in coordinate representation F^^ can be 
dropped. The disadvantage of this formula is that the FTs involve oscillating integrals 
which are numerically problematic. If the back-transformation F^^ is needed as in the 
case of Fo/5 and Fr it is better to perform the convolution directly. Similar to the FT the 
convolution can be reduced to the scalar case. The convolution of two scalar functions 
can further be reduced to a two dimensional integral: 

f*g{v) = ^j^fil~^^^_^^^J ^ dr f^dd f{r)g{^p^ (A.15) 

It is again evaluated with Gaussian integration methods. The second advantage is that 
there are no problems with slowly decaying functions. Sometimes there are large cancel- 
lations between different terms. In this case it is essential to use nonadaptive integration 
methods because they will not result in a loss of accuracy. 

The explicit reduction of the various correlators to the basic forms (|A.11|) and (|A.15|) 
is more or less trivial. The selfconsistency equation has been solved by iteration. The 
results are plotted in figure p.3| - |B.8| . 
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'"■A- 



£U4,X2l 



V-^Ov. 




^ 



Figure B.l: S is minimal at the river, small in the valley and large in the mountains, 
Therefore Z = Jdxidx2e~^^^^'^'^' is dominated by the valley. 
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Figure B.2: Constituent quark mass M{p) ~ P^^'^ip) in singular, regular and axial 
gauge for fixed instanton radius p in arbitrary normalization. For a given momentum the 
corresponding lowest curve may be interpreted as the "most physical" one. 
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Figure B.3: Pseudoscalar triplet correlator normalized to the free massless quark corre- 
lator. The pion coupling constant X^^ and the continuum threshold E^^ are fitted in order 
to match the spectral ansatz with the theoretical sum of the free and the connected part. 
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Figure B.4: Pseudoscalar singlet correlator normalized to the free massless quark cor- 
relator. There is a strong repulsion in this channel and no boundstate is formed. The 
theoretical curve is compared to a curve obtained from a pure continuum spectrum above 
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Figure B.5: Scalar triplet correlator normalized to the free massless quark correlator. 
There is a strong repulsion in this channel and no boundstate is formed. The theoretical 
curve is compared to a curve obtained from a pure continuum spectrum above Es- 
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Figure B.6: Scalar singlet correlator normalized to the free massless quark correlator. 
The a mass m„ and coupling X^ and the threshold E„ are obtained from a spectral fit. 
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Figure B.7: Axial vector correlator normalized to the free massless quark correlator. The 
triplet and singlet correlator are equal because the connected part has been neglected. The 
tti and /i mass, coupling and threshold are obtained from a spectral fit. 
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Figure B.8: Vector correlator normalized to the free massless quark correlator. The 
triplet and singlet correlator are equal because the connected part is zero. The p and u 
mass, coupling and threshold are obtained from a spectral fit. 
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